GLOBAL WELLPOSEDNESS IN THE ENERGY SPACE 
FOR THE MAXWELL-SCHRODINGER SYSTEM 
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Abstract. We prove that the Maxwell-Schrodinger system in R 3+1 
is globally well-posed in the energy space. The key element of the 
proof is to obtain a short time wave packet parametrix for the mag- 
netic Schrodinger equation, which leads to linear, bilinear and trilinear 
estimates. These, in turn, are extended to larger time scales via a boot- 
strap argument. 



1. Introduction 

The Maxwell-Schrodinger system in R 3+1 describes the evolution of a 
charged non-relativistic quantum mechanical particle interacting with the 
classical electro-magnetic field it generates. It has the form 

{iut — Aau = 4>u 
- Acft + d t divA = p, P=\u\ 2 
UA + V(d t (ft + divA) = J, J = 2Im(u, V A u) 

where u is the wave function of the particle, ((ft, A) is the electro-magnetic 
potential, 

(u,A,<ft) :l 3 xK-»Cxlx M 3 

and Va = V - iA, A A = V\. 

The system is invariant under the gauge transform: 

(u',cft',A') -► (e iX u,(ft-d t X,A + VX), X:R 3 xR^R 

where A:i? 3 xl^l. To remove this degree of freedom we need to fix the 
gauge. In this article we choose to work in the Coulomb gauge 

(2) divvl = 

Under this assumption, the system can be rewritten as: 

iut — A.au = <ftu 
OA = PJ 



(3) 



where (ft = (—A) 1 (|n| 2 ) and P = 1 — VdivA 1 is the projection on the 
divergence free vectors functions - also called Helmholtz projection. We 



The first author was partially supported by NSF grant DMS0738442. 
The second author was partially supported by NSF grant DMS0354539. 

1 



2 



IOAN BEJENARU AND DANIEL TATARU 



consider the above system with a set of initial data chosen in Sobolev spaces: 

(u(0),A(0),A t (0)) = (u ,A ,A 1 ) eH s xH°x H*' 1 

The gauge condition (2) is conserved in time provided the initial data 
(A ,Ai) satisfies it due to the form of the second equation in (3). 

The conserved quantities associated to the system are the charge and the 
energy, 

Q(u) = / |u| 2 dx 

E(u)= [ \V A u\ 2 + h\A t \ 2 + \V X A\ 2 ) + l\V<P\ 2 dx 
Jrs 2 2 

The local well-posedness of the system in various Sobolev spaces above 

the energy level is known, see [12], [9]. On the other hand the existence of 

weak energy solutions is established in [2]. The main outstanding problem 

which we seek to address is the well-posedness in the energy space. Our 

result is 

Theorem 1. The Maxwell-Schrddinger system (3) is globally well-posed in 
the energy space H 1 x H 1 x L 2 in the following sense: 

i) (regular solutions) For each initial data 

(u ,A),Ai) €H 2 xH 2 xH 1 
there exists an unique global solution 

(it, A) G C(R, H 2 ) x C(R, H 2 ) n C^R, H 1 ). 

ii) (rough solutions) For each initial data 

(«o,A),Ai) eH 1 xH 1 xL 2 

there exists a global solution 

(u,A) G CQR,^ 1 ) x C(R, H 1 ) n C 1 (R, L 2 ). 

which is the unique strong limit of the regular solutions in (i). 

Hi) (continuous dependence) The solutions (u, A) in (ii) depend continu- 
ously on the initial data in H 1 x H 1 x L 2 . 

We remark that in the process of proving the above results we establish 
some additional regularity properties for the energy solutions (u, A) which 
suffice both for the uniqueness and the continuous dependence results. Tra- 
ditionally these regularity properties are described using X s,b type spaces. 
Instead here we use the related U 2 and V 2 type spaces associated to both 
the wave equation and the magnetic Schrodinger equation. These are in- 
troduced in the next section; for more details we refer the reader to [7], [6], 
[3]- 

Remark 2. We note that in some directions our analysis yields stronger 
results than as stated in the theorem. Precisely, the same arguments as those 
in Section 5 also yield: 
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a) Local in time a-priori estimates in H 13 x H 1 x L 2 for (5 > \. This is 
exactly the range allowed for (3 in Lemma 25 and Lemma 26 (a). 

b) Local in time well-posedness in H 13 x H 1 x L 2 for (3 > |. This reduced 
range arises due to Lemma 26 (b). 

The nonlinearities on the right hand side of both equations in (3) are fairly- 
mild. Indeed, if A^ were replaced by A then it would be quite straightfor- 
ward to iteratively close the argument in X s,b or Strichartz spaces. For 
the magnetic potential A it is quite reasonable to hope to obtain an X s,b 
type regularity. Thus the main difficulty stems from the linear magnetic 
Schrodinger equation 

(4) iu t - A A u = f, u(0) = u Q 

The linear and bilinear estimates for L 2 solutions to (4) are summarized 
in Theorem 9 in Section 4. The rest of the section is devoted to the well- 
posedness of (4) in H 2 , H~ 2 and intermediate spaces. 

The proof of our main result is completed in the following section. The 
first step is to establish local in time a-priori bounds for solutions to (3), 
first in H 1 and then in more regular spaces. This is done by treating the 
nonlinearities on the right of the equations in a perturbative manner. The 
transition from local in time to global in time is straightforward, using the 
conserved energy. The second step is to establish the continuous dependence 
on the initial data. This is a consequence of a Lipschitz dependence result 
in a weaker topology. Precisely, we show that the corresponding linearized 
equation is well-posed in L 2 x Hz x H~2. 

The rest of the paper is devoted to the study of L 2 solutions for (4), with 
the aim of proving Theorem 9. Previous approaches establish Strichartz es- 
timates with a loss of derivatives for this equation in a perturbative manner, 
starting from the free Schrodinger equation. This no longer suffices for A 
in the energy space, and instead one needs to study directly the dispersive 
properties for the linear magnetic Schrodinger equation. Our approach uses 
some of the ideas described in [5] and [6]. 

To each dyadic frequency A, we associate the time scale A -1 . On this time 
scale we show that at frequency A the equation (4) is well approximated by 
its paradifferential truncation, which is roughly 

(5) iu t - & A< ^u = /, u(0) = u 

Following the ideas in [1], [8], in Sections 6, 7 we obtain a wave packet 
parametrix for the equation (5) on the A -1 time scale. This allows us to 
prove sharp Strichartz and square function estimates, as well as bilinear L 2 
bounds and trilinear estimates. 

In the last section we extend the linear, bilinear and trilinear estimates 
to larger time scales. A brute force summation of the short time bounds 
yields unacceptably large constants. Heuristically, the summation can be 
improved by taking advantage of the localized energy estimates for the mag- 
netic Schrodinger equation. However these are not straightforward. Our 
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idea is to obtain them from a weaker generalized wave packet decomposi- 
tion where the localization scales in position and frequency are relaxed as 
the time scale is iteratively increased. 

2. Definitions 

Throughout the paper we use the standard Lebesgue spaces L% and mixed 
space-time versions L\L% which are defined in the standard way. To measure 
regularity of functions at fixed time we use the standard Sobolev spaces H^,. 
Additional space-time structures will be defined in the next section. 

We now introduce dyadic multipliers and a Littlewood-Paley decomposi- 
tion in frequencies. Throughout the paper the letters A, fi, v and 7 will be 
used to denote dyadic values, i.e. A = 2 l for some i G M. 

We say that a function u is localized at frequency A if its Fourier transform 
is supported in the annulus {|£| G [§,8A]} if A > 2, respectively in the ball 
{|£|<8}if A = l. 

By S\ we denote a multiplier with smooth symbol s\{£) which is sup- 
ported in the annulus {|£| G [^,2A]} for A > 2 respectively in the ball 
{|£| < 2} if A = 1 and satisfies the bounds 

(6) \d a s x (0\ <c Q A-H 

By S < \ we denote a multiplier with smooth symbol s<a(£) which is sup- 
ported in the ball {|£| < 2A}, equals 1 in the ball {|£| < A/2} and satisfies 
(6). All implicit constants in the estimates involving S\, S < \ will depend 
on finitely many seminorms of its symbol, i.e. on c a for \a\ < N, for some 
large N. 

Associated to each A we also consider S\ to be a multiplier whose symbol 
s\ satisfies (6), is supported in {|£| G [f,4A]} and equals 1 in the support 
of { |£| 6 [^,2 A]}. The last condition implies that 

S\S\ = S\ 

Similarly we consider S\ to be a multiplier whose symbol satisfies (6), is 
supported in { |^| G [A/8, 8A]} and equals 1 in the support of {|£| G [A/4,4A]}. 

3. V 2 AND U 2 TYPE SPACES 

Let H be a Hilbert space. Let V 2 H be the space of right continuous H 
valued functions on R with bounded 2-variation: 

\\u\\y2 H = SUP V |K* i+ l) - U(U)\\ 2 H 

(ti)eT t 

where T is the set of finite increasing sequences in ML 
Let U 2 H be the atomic space defined by the atoms: 

i i 
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for some (tf) G T. We have the inclusion U 2 H C V 2 H but in effect these 
spaces are very close, and also close to the homogeneous Sobolev space Hz. 
Precisely, we can bracket them using homogeneous Besov spaces as follows: 

(7) b\ x CU 2 CV 2 C b| >00 

We denote by DU 2 H the space of (distributional) derivatives of U 2 H 
functions. Then there is also a duality relation between V 2 H and U 2 H, 
namely 

(8) (DU 2 H)* = V 2 H 

For more details on the U 2 and V 2 spaces we refer the reader to [7] and [3] . 
Given an abstract evolution in H, 

iut = B(t)u, u(0) = uq 

which generates a family of bounded evolution operators 

S(t,s):H^H, t,seR 

we can define the associated spaces U^H, V^H and DU^H by 

\H\viH = \\S(0,t)u(t)\\ V 2 H , \\u\\ uljH = \\S(0,t)u(t)\\ v2H 

respectively 

DUlL 2 = {(id t - B)u; u G U 2 B H] 

On occasion we need to compare the above spaces associated to closely 
related operators. For this we use the following 

Lemma 3. Let H be a Hilbert space and B(t), C(t) two families of bounded 
selfadjoint operators in H . Suppose that for e sufficiently small we have 

- C)u\\ DU 2 H < eH^h 

Then 

\\ u \\u%H ~ IMIvjtf ~ ll^lly^i?! \\f\\DU%H ~ UWdu^h 



Proof. By conjugating with respect the B flow we can assume without any 
restriction in generality that B = 0. Then we can solve the equation 

iut — Cu = 0, u(0) = uq 

by treating C perturbatively to obtain a solution 

u = u + eu 1 (t), \\ui\\u2 H < \\u \\ H 

Applying this to each step in Uq atoms we obtain 

\H\lPH ^ \\ u \\u*H 

for arbitrary u. 
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For the converse, applying the above result to each step in a U 2 H atom 
we conclude that for each u G U 2 H we can find u\ G U 2 H so that 

\\u + euiWfjijj + ||«i||[72 H < \\u\\u2 H 

Iterating this shows that 

IMIt/gH ^ IMIt/ 2 .ff 

Hence U%H = U 2 H. 

Consider now / G DU^H. Then / = iu t - Cu for some u G = U 2 H. 

Since C maps ?7 2 i7 to DU 2 H this implies that / G DU 2 H. Conversely, if 
/ G DU 2 H then we can solve the inhomogeneous equation 

iu t -Cu = f, u(0) = 

iteratively in U 2 H. This gives a solution tx G C/ 2 // = U^H, therefore 
/ G DU^H. We have proved that DU 2 H = DU^H. Then the last relation 
V 2 H = V 2 H follows by duality. □ 

Following the above procedure we can associate similar spaces to the 
Schrodinger flow by pulling back functions to time along the flow, namely 

II II II it A ii ii ii ii it A ii 

IMIvj;L 2 = ll e w llv 2 L 2 > IMIt/2 L 2 = ||e U\\ V 2 L 2 

The magnetic Schrodinger equation has time dependent coefficients, so we 
replace the above exponential with the corresponding evolution operators. 
We denote by S (t,s) the family of evolution operators corresponding to 
the equation (4). These are L 2 isometries. Then we define 

IMIV|L2 = \\S(0,t)u(t)\\y2 L 2, |M|[/2 L 2 = \\S(0,t)u(t)\\ V 2 L 2 

These spaces turn out to be a good replacement for the X ' 2 space associated 
to the Schrodinger equations. We also define 

DU 2 A L 2 = {{id t - A A )u; u G U 2 A L 2 } 

By (8) we have the duality relation 

{DU 2 A L 2 )* = V\L 2 

When solving the equation (4) we let / G DU\L 2 , and we have the 
straightforward bound 

(9) IMIc/ 2 L 2 < INIIl 2 + II/IIdc/ 2 l 2 

In our study of nonlinear equations later on we need to estimate multilinear 
expressions in DU\L 2 . By duality, this is always turned into multilinear 
estimates involving V^L 2 functions. 

Finally, we define similar spaces associated to the wave equation. The 
wave equation is second order in time therefore we use a half-wave decom- 
position and set 

||„, II _ M_iitl.DL.il ||,, || _ ||„iit|DL.|| 

IMIv 2 L 2 — ll e u \\V 2 L 2 i IMIl/ 2 L 2 — ll e U \\V 2 L 2 
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Then the spaces for the full wave equation are 

IMIt/^L 2 = IMI[/2 L 2 +l /2 L 2, lkllv^L2 = |M|y2 L 2 + y2 i2 

For the inhomogeneous term in the wave equation we use the space DU^/L 2 
with norm 



DU^L 2 - \\J\\DUlL 2 nDUlL 2 

Then to solve the inhomogeneous wave equation we use 

ll Vw lli/^L2 ^ l|Vu(0)|| L 2 + IpuWuu^^ 

Similarly we set 

IMIt/2^ = ll(^) S «||(72L2+t/2 L 2, IMIv^H" = \\( D x) S u\\ V 2 L 2 +V 2 L 2 

and 

\\f\\DU* w Hs = \\(D x ) S f \\ DU 2 L 2 

Such spaces originate in unpublished work of the second author on the wave- 
map equation, and have been successfully used in various contexts so far, 
see [7], [6],[1],[3]. The Strichartz estimates for the wave equation turn into 
embeddings for U^rH s spaces. If the indices (p, q) satisfy 

(10) - + - = l, 2<p<oo 

p q 

then we have 

(11) IMUpl* < IMi^j 

If we consider frequency localized solutions to the wave equation on a 
very small time scale, then the wave equation is ineffective. Precisely, 

Lemma 4. Let B < \ be a function which is localized at frequency < A. Then 

(12) ll^<A||;72 r (/;L2) W ||£<a||c/2(/;L2), |/| < A _1 

Proof. This follows from the similar bound for the corresponding half-wave 
spaces U±L 2 , and by Lemma 3 it is a consequence of the fact that for a 
short time the spatial derivatives in the half wave equation can be treated 
perturbatively, 

IIIAz|£<aIIdi/2(7,l2) < HIAeI^aIIli^l 2 ) ^ ^II- b <aIIl i (/,l 2 ) 

< II^aIIl^^l 2 ) ^ |^|A ||-B<a||c/2(/,l2) 

□ 

The finite speed of propagation for the wave equation allows us to spatially 
localize functions in the U^L 2 spaces. For R > we consider a covering 
{Qf)ieI A °f ^ 3 with cubes of size R. Let xf be an associated smooth 
partition of unity. Then we have the following result: 
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Lemma 5. Assume that I is a time interval with \I\ < R. Then: 

Proof. By rescaling we can take R = 1. Without any restriction in generality 
we can also assume that |/| =1. We prove that the result holds for one of the 
two half-wave spaces, say U+(I; L 2 ). It is enough to verify (13) for atoms, 
and further for each step in an atom. Hence we can assume that u solves 
the half wave equation 

(id t + \D x \)u = 

Then we have 

(id t + \D x \)(xfu) = [\D x \, X f]u 
By standard commutator estimates we have at fixed time 

Y,\\[\D x lxf]u(t)\\h<\Ht)\\h 

Then 

E WxM 2 ulriI;L2) < E llxf^(0)||i 2 + \\[\D x \,xf]u{t)\\W, < ||n(0)||| 2 

□ 

Due to the atomic structure, in many estimates it is convenient to work 
with the U 2 type spaces instead of V 2 . In order to transfer the estimates 
from U 2 to V 2 we use the following result from [3]: 

Proposition 6. Let 2 < p < oo. If u £ V 2 H then for each < e < 1 there 
exist u\ G U 2 H and 112 £ U P H such that u = u\ + U2 and 

(I 4 ) I Iner^killc/ 2 ^ + £~ 1 \\u2\\uph < IMIv^H 

Here U p is defined in the same manner as U 2 but with the I 2 summation 
replaced by an V summation. One way we use this result is as follows: 

Corollary 7. Let e > and N arbitrarily large. Then 

V^H S C U^H s - e + L°°H N 



Following is another example of how this result can be applied. Typically 
in our analysis we prove dyadic trilinear estimates of the form 
r3 



(15) 



/ / Sx^SxzV, S\ s Bdxdt 
J 1 Jr 



< Cl(l^|,Al23)||^||[72L2||«||[/2 L 2||B|| [/ 2 r2 



What we need instead is an estimate with one U 2 replaced by a V 2 , say 



(16) 



J 1 Jr 



uS\ 2 v, S\ 3 Bdxdt 



< C 2 {\I\,\i2?,)\\u\\ U 2 aL 2 \\v\\jj2 aL 2 \\B\\ V ^ L 2 
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Denoting 

Amax = max{Ai, A2, A3}, 
due to Proposition 6 we can easily show that 

Lemma 8. Assume (15) holds and \I\ < 1. Then (16) holds with 

C2(\I\, A123) = Cl(|^|) A123) hi Xmax- 

The same holds if the V 2 structure is placed on any of the other two factors 
in (16). 

Proof. Without any restriction in generality we assume that Ai, A2 and A3 
are so that the integral in (15) is nontrivial. Taking u, v and B to be time 
independent frequency localized bump functions we easily see that 

(17) C-i(|J|,Ai 23 )>|I|A^ 

for some sufficiently large N . 

For each < e < 1 we decompose B = B\ + B2 as in Proposition 6. For 
B\ we use (15) while for B2 we use Bernstein's inequality to estimate 
r-3 

~ I^IAmaJMlL^L^MlL^LHI^HL^L 2 



/ / S\ 1 uS\ 2 v,S\ 3 B 2 dxdt 
J 1 Jr 



^ l^|A ma a; IMI[/ 2 L 2 \\ v \\u\ L 2 1 1 ^2 1 1 L 2 



Adding the B\ and the B2 bounds gives 

CMI4Ai23) <|lne|C 1 (|/|,A 1 23) + e |/|Al a; 
We set e = A~^. Then the conclusion of the Lemma follows due to (17). 

□ 

4. The linear magnetic Schrodinger equation 

In this section we summarize the key properties of solutions to the homo- 
geneous and inhomogeneous linear magnetic Schrodinger equation 

(18) iu t - A A u = 0, u(0) = u 

(19) iut - A A u = f, u(0) = u 

in L 2 , and we use them in order to show that the above equation is also 
well-posed in H 2 , H~ 2 and in intermediate spaces. 

We assume that A G U^H 1 with = 0. All constants in the estimates 
depend on the U^H 1 norm of A which is why we introduce the notation 
X < A Y, which means X < C(\\A\\ u2 H i)Y. 

The trilinear estimates are concerned with integrals of the form 

rT 



/ Ai,A2,A 3 ( u >' y > B ) = / / S Xl uSx 2 v S Xs B dxdt 
Jo Jr 3 



where u and v are associated to the magnetic Schrodinger equation and B 
is associated to the wave equation. In order for the above integral to be 
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nontrivial the two highest frequencies need to be comparable. Thus by a 
slight abuse of notation in the sequel we restrict ourselves to the case 

{Ai, A 2 , A 3 } = {A,A,/4, M<A 
With these notations we have 

Theorem 9. For each A £ U^H 1 with V • A = the equation (18) is 
well-posed in L 2 . For each e > there exists 5 > so that the following 
properties hold for < T < 1 : 

(i) Strichartz estimates: 

£ i 2 3 3 

(20) \\S X U\\ L P(0,T;L1) <A Tv \p \\u\\jj2 L 2 - + - = -, 2 < p < oo 

p q z 

(ii) Local energy estimates. For any spatial cube Q of size 1 we have 

(21) \\Sxu\\ L 2^ T . Q) < A T 5 \-^\\ u \\ u2aL2 . 

(Hi) Local Strichartz estimates. For any spatial cube Q of size 1 we have: 

( 22 ) \\ S \ u h 2 (0,T;L6(Q)) <A T S X e \\u\\ u ^ L 2. 

(iv) Trilinear estimates. For any < T < 1 and /i < A we have 

(23) \llxju, v, B)\ <4 T 5 X e min (1, fiX^M^ ||S||^ i2 . 
On the other hand if fi <^ A then 

(24) \ll^ x (u,v,B)\ ^a^A-H^IIuII^II^II^IIBII^ 

The proof of this theorem is quite involved and is relegated to Sections 6- 
10. 

The smallness given by the T s factor is needed in several proofs which 
use either the contraction principle or bootstrap arguments. However, this 
factor is nontrivial only in (20). Indeed, we have 

Remark 10. Assume that the conclusion of Theorem 9 holds without the 
T 5 factor in (21), (22), (23) and (24). Then the conclusion of Theorem 9 
holds in full. 

Proof. For (21) we observe that 

\\ S \ u \\l*(0,T;Q) ^ T ^IMIl°°L 2 ^ \\u\\ulL 2 - 

Interpolating this with (21) without the T s factor yields (21) with a T 5 
factor. By Bernstein's inequality the same argument works for (22). 
For (23) we can also write the obvious estimate 

I^K"' 5 )! ~A T^5|| u || LooL2 || v || Ltx)L2 ||5|| L ^ L 2. 

which is then interpolated with (23) without the T s factor. The same argu- 
ment applies for (24). □ 



MAXWELL-SCHRODINGER 11 

Next we turn our attention to the H 2 and H~ 2 well-posedness for (18). 
We make the transition from L 2 to H 2 and H~ 2 using the coercive elliptic 
operator 1 — A^. Its properties are summarized in the following 

Lemma 11. For each < s < 2 the operator 1 — is a diffeomorphism 

1 - A A : H s -> H s ~ 2 

which depends continuously on A £ H 1 . 

The proof uses standard elliptic arguments and is left for the reader. 

Using the above operator we define the spaces U\H 2 , V\FL 2 , respectively 
DU\H 2 by 

\\ U \\UIH* = IK 1 - A A)W|| 1 /2 L 2, \\u\\y2 H 2 = ||(1 - A A ) U \\ V 2 L 2 

respectively 

ll/llw^ = |l(l-AA)/|UiL 2 

Remark 12. The reason we use the U , V notation above is to differentiate 
these spaces from the U\H 2 , V\H 2 , DU\H 2 spaces which should be defined 
as in the previous section, with respect to the H 2 flow of (18). This is not 
possible at this point, as we have not yet proved that (18) is well-posed in 
H 2 . However, after we do so we will prove that the above two sets of norms 
are equivalent. After that the U, V notation is dropped. 

We can transfer the estimates from Theorem 9 to the U\H 2 spaces by 
making an elliptic transition between A^ and A: 

Lemma 13. Let A G U^H 1 with V • A = 0. Then for each e > there 
exists 5 > so that the following properties hold: 

(i) Forp,q as in (20) we have the Strichartz estimate 

- -2+- 

(25) \\S\u\\ LP{0jT . Lq) < a Tp\ p||u||^2 H 2 

(ii) Elliptic representation. Each u £ U\H 2 can be expressed as 
(26) 

U=(l-A) 1 (u e + U r ), \\Ue\\ulL^ + T 6 \\ u A\l^{0,T;H^) Si 

(Hi) Local energy estimates. For any spatial cube Q of size 1 we have: 

(27) \\Sxu\\ L 2 ([0>T]xQ) <4 T 5 \- 2 -^\\u q \\ 2 aH 2 

(iv) Local Strichartz estimates. For any spatial cube Q of size 1 we have: 

(28) ||5'a'"||l2( 0)T . L 6(q)) <AT 5 \~ 2+e \\u G \\fj2 AH 2 
Proof, (i) The Strichartz estimate (25) follows from 

(29) ||5 A (1 - A)u\\ L P Lq < A TpAp||u||^ 2 
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We use the identity 

(30) (1 - A)u = (1 - A A )u - 2iAX7u - A 2 u = (1 - A A )u - R A {u) 

and estimate each of the three terms. From the definition of U\H 2 and (20) 
we have: 

||Sa(1 - A A )u\\ LPLq < A Tv\p ||u||y-2 H 2 

For the second term we use Bernstein's inequality and the exponent relation 
in (20) to estimate 

\\S x {AVu)\\ LPLq <4 Thhs x {AVu)\\ L oo L 2 

<Tp\p\\A\7u\\ i 
l i. 

< T? \p \\A\\ L oo H i\\u\\ L oo H 2 
Similarly for the last term we obtain 

\\S x {A 2 u)\\ LPLq < TpXp\\A\\ 2 looh1 \\u\\ L oo H 2 

This concludes the proof for (29) which implies (25). 
(ii) By (30) we can set 

u e = (1 — A A )u, u r = —R A (u) = —2iAX7u — A 2 u 

Hence it remains to prove that 

(31) ||S A (AV«)|| L2 + \\S x (A 2 u)\\ L 2 <4 T 5 A- 1+e ||n||^^ 

By the argument in Remark 10, here and for the rest of the proof of the 
lemma we can neglect the T s factors. 

We decompose the expression S\(AVu) as 

(32) Sx(AVu) = S x I s i As ^ u + ^ SxAS^Vu + ^ S 7 AS 7 Vi 

\7<A 7<A 7>A 

For exponents (p,q) satisfying (10) we invoke the Strichartz estimates (11) 
for the wave equation. By Bernstein's inequality and (25) we can derive a 
similar bound for u, 

_ 2 

(33) II^ViiHi/iip < A 7 p \\u\\u^ H 2- 

The L 2 bound for the product is obtained by multiplying the last two in- 
equalities. For the first term in (32) we take q close to oo, for the second we 
take p = oo, while for the third any choice will do. 

For the expression Sx(A 2 u) we take a triple Littlewood-Paley decompo- 
sition, 

S x {A 2 u) = ^ S x (A Xl Ax 2 ux 3 ) 

Ai,A2,A3 

Then we must have X max > A. 
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If ^max = A3 then we use the above Strichartz inequalities and Bernstein's 
inequality to estimate the triple product as 

||4\i A\2 U A 3 1 1 £2 < \\A Xl ||l4l°o \\A\ 2 Hl 4 L°° 1 1 "A3 IIl°°L 2 
1 1 

<a \? A 2 4 A 3 2 || A Xl Wu^m Pa 2 Wu^hi IK 3 \\u\H* 

The summation with respect to Ai, A2, A3 is straightforward. 

If A3 <C X max then there are two possibilities. One is X max = A, in which 
case we assume w.a.r.g. that A = Ai > A2,A3 and estimate 

ll^-Ai ^A 2 W A 3 Hi 2 < II^AllLPL^I^Aalli^ll^Aallz 00 

_ 1 —'A _ 1 

< A X 1 P X 2 9 A3 ^I^Aillf/^ifill^Azllt/^ifillwAslli/Z^ 

with q close to infinity. The other possibility is X max » A, in which case we 
must have Ai = A2 3> A3. Then we estimate the triple product as above, but 
the choice of p and q is no longer important. The proof of (31) is concluded. 

(iii) We use the representation in (26). The bound for u r holds without 
any localization. For u e we can write 

S x (l - A)-V = A~ 2 (A 2 (1 - A^SJSxUe 

where the symbol of S\ is still supported at frequency A but equals 1 in 
the support of the symbol of S\. The operator X 2 A~ 1 S\ is a unit mollifier 
acting on the A -1 scale, therefore it is bounded in IqL 2 ([0, 1] x Q). 

(iv) We use the representation in (26). By Bernstein's inequality the 
bound for u r holds without any localization. For u e we argue as above. □ 

Next we define similar spaces U A H~ 2 , V A H^ 2 , respectively DU\H~ 2 in 
a manner similar to the H 2 case, namely 

IMIj/2^-2 = IK 1 - A ^r lu llt/2 L 2, M|y2ff-2 = 11(1 - A A ) _1 u||y2 L 2, 

respectively 

Wf\\DU*H-* = IK 1 " ^A)~ l f\\ D uiL* 
Due to the duality relation 8 we have the H 2 — H~ 2 duality 

(34) {DU 2 A H- 2 )* = VlH 2 , (DU 2 A H 2 )* = V\R- 2 

For functions in U\H~ 2 we can prove results similar to Lemma 13: 

Lemma 14. Let A G U^H 1 with V • A = 0. Then for each e > there 
exists 5 > so that the following properties hold: 

(i) For p,q as in (20) we have the Strichartz estimate 

- 2+- 

(35) \\S\u\\ LP ^ T . Lq) < A TpX r\\u\\jj2 H - 2 

(ii) Elliptic representation. Each u G U\H 2 can be expressed as 
(36) 

U=(l- A)(u e + U r ), \\Ue\\u\L2 + T S \\ u r\\L2(0,T;m-*) <A \\u\\u2 H -2 
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(Hi) Local energy estimates. For any spatial cube Q of size 1 we have: 

(37) \\S\ u \\l 2 ([o,t]xQ) ~a T 5 A 2- 2 +€ ||uo||j/-2 H -2 

(iv) Local Strichartz estimates. For any spatial cube Q of size 1 we have: 

(38) \\S\u\\ L 2 (0)T . L6m < A T s X 2+e \\u \\(j^ H -2 

Proof. The proof is similar to the proof of Lemma 13, so we merely outline 
it. We denote v = (1 — Aa)~ 1 u. Then 

u= (1 - A)v- R A (v) 

To prove (35) we use (20) for (1 — A)v and it remains to show that 

\\SxV(Av)\\ LPLq + \\S x (A 2 v)\\ LPLq <4 T 5 \ 2+1 v\\v\\ ulL 2 

which is obtained using only the energy estimates for A and v. 
For (36) we set 

u e = v, u r = (1 — A) -1 ^^) 
Then it remains to show that 

\\S x V(Av)\\ L 2 + \\S x (A 2 v)\\ L . <4 T^IMI^ 

But this is obtained in the same manner as (31) from the Strichartz estimates 
for A and U2- 

Finally, (37) and (38) are proved exactly as in the previous lemma. □ 

Next we turn our attention to the trilinear bounds, namely the H 2 and 
H~ 2 counterparts of (23) and (24). For uniformity in notations we set 
U 2 A L 2 = U\L 2 . Then 

Lemma 15. Let k,l G {—2,0,2}. Then for any < T < 1 and /j, < X we 

have 

(39) |/£a,>, v, B)\< A T s X e min (1, fiX'^X^-i^^ | M | ||£||^ L2 
while if fi <^ X then 

(40) \Ll^ x (u,v,B)\ ^A^A-H^^-'A-^lull^llvll^HBll^ 

Proof. The T s factor can be neglected by the argument in Remark 10. We 
represent 

k 

u = (1 — A)~ 2 (u e + u r ) 

where u e , u r are chosen as in (26) if k = 2, as in (36) if k = — 2 and with 
u r = if k = 0. Similarly we set 

v = (1- A)-i(v e + iv) 

We begin with (39) and consider all four combinations. The estimate for 
u e and v e is exactly (23). The estimate for u e and v r reads 

I^A^K, V r ,B)\ < A X e min (1, /iA - 5)||«e||l^La KIIl 2 ^- ll^llt/^ 
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and is a consequence of the stronger bilinear L 2 estimate 

(41) \\S X U e S^B\\ L 2 < A ^X^lUeWjj^ ||-B||[7^ L 2 

Due to the finite speed of propagation for the wave equation, see (13), we 
can localize this to the unit spatial scale. But on a unit cube Q we use the 
local Strichartz estimate (22) for u e and the energy estimate for B. 

The estimate for u r and v e is similar. The estimate for u r and v r reads 

\ I l,\A Ur > Vr > B }\ ~ A ^ emin ( 1 >^~^)IKIIl^Hi-«l|Wr||LaHi-«||-B|| £ ^ L 2 

and is easily proved using L 2 bounds for S\u r , S\v r and an L°° bound for 
S^B. 

The proof of (40) is similar. For later use we note the bilinear L 2 esti- 
mates, namely 

(42) WSpVeSxBWv < A ^ +e Kb»L2 \\ B \\u^ 
respectively 

(43) \\S^S x u e S x B)\\ L 2 < A A^KII^II^IIt/^ 

Both are proved by localizing to a unit spatial scale and then by combining 
the local Strichartz estimate (22) for u e and v e and the energy estimate (11) 
for B. 

□ 

Now we consider the H 2 well-posedness of (18). 

Proposition 16. Let A G Uyy with V ■ A = 0. Then the equation (19) is 
well-posed in H 2 . In addition we have 

U 2 A H 2 = U 2 A H 2 , VlH 2 = VlH 2 , DU\H 2 = U 2 A H 2 

with equivalent norms. 

Proof. In order to solve the equation (18) with initial data uo G H 2 we 
consider the equation for v = (1 — Aa)u which has the form 

(id t - A A )v = 2(A t V - iAA t )u 

Expressing u in terms of v we obtain 

(44) (id t - A A )v = 2(A t V - iAA t )(l - Aa^v 

We seek to solve this equation perturbatively in U\L 2 . For this we need 
first to establish suitable mapping properties for the operator AtV — iAAt. 

Lemma 17. The operator A t V — iAA t satisfies the space-time bound 

(45) \\(A t V - 2iAA t )u\\ DUlL2 < A T 5 ||u||^ 2 
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Proof. By duality, (45) follows from the bounds 



(46) 
respectively 
(47) 



BVuvdxdt 



o JR 3 



<aT* 



\ u \\u*H*W v \\v%L 2 \\ B \\u% r L i 



fj 




Jo Jl 


R 3 



<a T d \\B\ 



u&l*\\Mu&h4 u \\u*h*\\ v \\v*l* 



fj 


BVuvdxdt 


Jo Jl 


R 3 



W •'"w** " " ^ A 1 

To prove (46) we use a triple Littlewood-Paley decomposition to write 

< E v ' b ) + E 7 J a,a(v«, 5) 

/i<A M<A 
+ E J A T M,A(V^^i?) 

Then for each term we use the corresponding bounds (39), and (40) with a 
In A correction coming from the use of Proposition 8. The summation with 
respect to [i and A is straightforward. 

For (47), by (13) the norms of A and B are I 2 summable with respect 
to unit spatial cubes. Hence without any restriction in generality we can 
assume that both A and B are supported in a unit cube Q. For A we use a 
Strichartz estimate, for B the energy and for u a pointwise bound. Finally, 
for v we interpolate the local energy estimates with the local Strichartz 
estimates to obtain 



\S\v\\ „ 12 <a A 4~ 
1 A u l 2 l~f ~ 



which leads to 



Then we can estimate 

rT f * 

/ ABuvdxdt 

'0 JQ 



Ml „ 12 S,a IMK/2/-2 

1 n L 2 LS- ~ A 11 nV A L 



< r4||S|| LO o L 2||A|| £ 4 £ 12||«||l,oo||t;|| i2i ^ 



<a Ti\\B\ 



□ 



We now return to the equation (44). By the definition of the U\L 2 norm 
and (45) we have 



(48) 



||2(A t V — iAA t ){l — Aa) 



By (9) it follows that we can solve (44) perturbatively in U\L 2 on short 
time intervals. This gives a solution u = (1 - Aa^v £ U 2 A H 2 for (19). 
Furthermore, we obtain the bound 

Kt)-s(t,oMo)|| £2 < A T 5 lk(o)|| £2 
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where S(t, s) is the evolution associated to (18). In particular this shows that 
v G C([0, 1]; L 2 ). An elliptic argument allows us to return to u and conclude 
that u G C([0, 1]; H 2 ). This concludes the proof of the H 2 well-posedness. 

Finally, we show that U A H 2 = U\H 2 and the other two similar identities. 
Via the operator / — A A these two spaces can be identified with the U 2 L 2 
spaces associated to the equations (18), respectively (44). But by Lemma 3, 
these are equivalent due to (48). 

□ 

Next we consider the well-posedness in H~ 2 , which is essentially dual to 
the H 2 well-posedness. 

Proposition 18. Let A G U^H 1 with V • A = 0. Then the equation (19) 
is well-posed in H~ 2 . In addition we have 

U 2 A H- 2 = U A H' 2 , V\H~ 2 = VlH~ 2 , DU 2 A H- 2 = U 2 A H- 2 

with equivalent norms. 

Proof. By Lemma 11 we can write the initial data uq as 

n = (1 - A A )v , v G L 2 

Then we seek the solution u for (18) of the form u = (1 — A A )v. The 
equation for v is 

(49) (idt - A A )v = 2(1 - A A )-\A t V - iAA t )v 
To solve it we need the following counterpart to (45): 

Lemma 19. The operator A t V — iAA t satisfies the space-time bound 

(50) \\{A t V - 2iAA t )u\\ DUlH - 2 < A T 5 ||u||^ L2 

Proof. By duality, (50) follows from the bounds 



T 



(51) 



respectively 
(52) 



BVuvdxdt 

R3 



<A T \\uWjj2 L 2\\V\\ V 2 H 2 \\BWjj2 t2 



ABuvdxdt 



T^lBWu^L^WMu^m \\ u \\u\ lz\\ v II vim 



These are almost identical to (46) and (47), and their proofs are essentially 
the same. □ 

The bound (50) allows us to solve the equation (49) perturbatively in 
U A L 2 and obtain a solution v G C([0, 1];L 2 ). This implies the H~ 2 solv- 
ability for (18). The second part of the proposition follows again from 
Lemma 3. □ 
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Having the well-posedness result in H 2 and H~ 2 allows us to prove well- 
posedness in a range of intermediate spaces. Given a positive sequence 
{m(A)} A=23 satisfying 

(53) 0<c<^-<C 

m(A) 

we define the Sobolev type space H(m) with norm 

\M\ 2 H{m) =^m\X)\\Sxu\\ 2 L2 
x 

The standard Sobolev spaces H a are obtained by taking m{\) = X a . 

We consider the solvability for (18) in H(m) under a stronger condition 
for m, namely 

(54) i<=«<4 

4 m(A) 

This guarantees that H(m) is an intermediate space between H~ 2 and H 2 . 

We need to describe H(m) in terms of H~ 2 and H 2 . To measure functions 
which are localized at some frequency A we can use the norm 

r\ /lii iiO /In iiO 

\\ u \\h x = ^~ IrIIj/z + A ||u|| h -2 

Ideally we would like to represent H(m) as an almost orthogonal superpo- 
sition of the H\ spaces with the weights m(\). However, this does not work 
so well if H{m) is "close" to either H~ 2 or H 2 . Instead we need to select a 
subset of dyadic frequencies which achieves the desired result. We denote 

moo = lim A 2 m(A) 

On 2 N U {oo} we introduce the relation "-<" by 

A -< fi 44> 2m(p) > m{\)(\ 2 ^r 2 + /U 2 A~ 2 ), fi, A < oo 
respectively 

oo -< ^ 44> 2m(/i) > [i^rnoo, fi < oo 

Definition 20. We say that a subset A(m) C 2 N U {oo} is m-representative 
if (i) for each G 2 N there exists A G A(m) so that A -< and (ii) for each 
fi G A(m) there is at most one A G A(m) \ {fi} such that A -< /i. 

Lemma 21. // m satisfies (54) i/ien an m-representative set A(m) exists. 
In addition, for each fi 6 2 M and K £ N we have 

(55) |{A G A(m); 2 K m(/x) > m(A)(A 2 ^ 2 + ,u 2 A~ 2 )}| < 4(if + 4) 
Proof. For each A G 2 N U {oo} we denote 

/ A = {/.G2 N U{oo}, A^} 
Due to (54) it is easy to see that I\ is an interval, 

/ A = [A-,A+], A~<A<A+ 
and the endpoints A~ and A + are nondecreasing functions of A. 
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We construct the set A(m) as an increasing sequence {Xj} in an iterative 
manner. Ao is chosen maximal so that Ao -< 1. Iteratively, Aj+i is chosen 
maximal so that Aj+i -< 2A^~. Either this process continues for an infinite 
number of steps, or it stops at some step k with A& = oo. The former occurs 
if moo = oo and the latter if vn^ < oo. 

The property (i) in Definition 20 is satisfied by construction. For (ii) 
we observe that Aj+i > 2A^ therefore Xj -ft Aj+i. On the other hand by 
construction we have Aj+2 7^ Xj. 

For (55) suppose fj, < Xj. Then 

m(ji) > -m(Xj)n 2 Xj 2 > -m(X j+2 )n 2 Xj 2 2 >■■■> 2" x " 2 m(A j+2 x)^ 2 Xj 2 2K 

A similar bound holds if we descend from fj,, and the conclusion follows. □ 

Since we allow 00 £ A(m) we need the equivalent of H\ in that case, 
which is defined by = H~ 2 . We note that at the other extreme we have 
H ± = H 2 . 

Lemma 22. Let m satisfy (54), and A(m) be an m-representative subset of 
2 JV U{oo}. Then 

\\ u \\H(m) ~ inf { Yl m ( A ) 2 H M A||H A , U= Yl U ^ 
AeA(m) AeA(ra) 

Proof. By Definition 20 we have a finite covering of 2 N with intervals 

2 N C (J h 

AeA(m) 

We consider an associated partition of unity in the Fourier space, 

AeA(m) 

For ji € I\ we have m{ji) ps m(X)(/j, 2 X~ 2 + X 2 fi~ 2 ) therefore we obtain 
||XA(ArHU( m ) ~ m(A)||x A (£ > x)«|| J H' A 

and the ">" inequality follows. For the reverse we use (55), which shows 
that the series Yl u >> ls almost orthogonal in H(m), 

{ux^Ux^Him) < 2 H *~ J U(A^m(A^|KJif A .|KJ.ff A . 

□ 

Finally we consider the well-posedness of (18) in H(m). 

Proposition 23. a) Assume that the sequence m satisfies (54). Then the 
equation (18) is well-posed in H{m). 

b) Furthermore, for each u G U\H(m) there is a representation 

u= £ Ux 

AeA(m) 
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with 

(56) Yl m '( A ) ( A " 4 H u a||^ 2 + A 4 ||u A ||^_ 2 ) < A \\u\\lj, H{m) 

AeA(m) 

c) The following duality relation holds: 

(DU\H(m))* = VlH{m- 1 ) 

Proof, a) We consider a dyadic decomposition of the initial data 

uo= Xx(D x )u 

AeA(m) 

and denote by u\ the solutions to (18) with initial data S\Uq. Then 

u = ^ U\ 

AeA(m) 

We can measure u\ in both H 2 and H~ 2 , 

^ 2 \\ u x\\c([0,l],H' 2 ) + ^ 2 \\ u x\\c([0,l],H-' 2 ) <A II^A^ollia 

After summation this gives 

Y m 2 (A) {\~*\\u x \\c(io,i],H*) + A 4 ||«a||c([o,i],//- 2 )) \\uo\\ 2 H (m) 

AeA(m) 

which, by (54), implies that u £ C([0, 1], H{mj) and 

IHIc([0,l],ff(m)) <A ||^o||ff(m) 

b) It suffices to consider the case when u is an U\H(m) atom. Then 
we consider a decomposition as in part (a) for each of the steps, and the 
conclusion follows. 

c) This is a direct consequence of (8). □ 

Finally, we can interpolate the properties from U\H 2 and U\H~ 2 to 
obtain properties for U 2 H(m). Indeed, we have the following 

Proposition 24. Let A G U^H 1 with V • A = 0. Assume that m, mi, mi 
satisfy (54) . Then the following estimates hold: 

(i) Forp,q as in (20) we have the Strichartz estimate 

- -l - 

( 57 ) \\S\u\\ LP ( 0tT . Lq) < A Tpm(\) X^WuWu^h^ 

(ii) Local energy estimates. For any spatial cube Q of size 1 we have: 

(58) \\Sxu\\ L 2 {[0 , T]xQ) < A T <5 m(A)~ 1 A~5 +e ||ti ||^2 /i -( m - ) 

(Hi) Local Strichartz estimates. For any spatial cube Q of size 1 we have: 

(59) \\Sxu\\ LH0:1 . L e m <4 T 5 m(\y 1 X e \\ 
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(iv) Trilinear estimates. For fj, < A we have 

tT r>w^ T s X e min (1, 2 ) , 

(60) \I XXll (u,v,B)\< A mi(A)m2(A) \Hu%H( mi )\MulH( m2 )\\B\\u^ 

w/u/e if jj, <C A i/ien 

(61) |/^, A (n, W)J B)| < A mi(A)m2(/u) II^II^H( mi )^II^H( m2 )PII^L 2 

Due to the representation in Proposition 23(b) this result is a straightfor- 
ward consequence of the similar results in H 2 and H~ 2 . 

5. Proof of Theorem 1 

We first establish an a-priori estimate for regular (H 2 ) solutions of (3). 
For this we need to consider the two nonlinear expressions on the right hand 
side of (3). We begin with the Schrodinger nonlinearity: 

Lemma 25. For j3 > \ and m satisfying (54) we have 

(62) \\M\DUl(0,T;H(m)) <A T S \\v || ^2 H{m) \\ u\\ ^2 m 

Proof. By duality the above bound is equivalent to the quadrilinear estimate 

nA"Vl U 2 )u 3 U4dxdt Sl^lKllt/^IKIIf^ \H0\\ u 3\\u 2 A H(m)\\ u 4\\ V 2 H( ±) 
,3 A A A A \ m > 

After a simultaneous Littlewood-Paley decomposition of the three factors 
A _1 (uiU2), U3 and 114 we need to consider the following three sums: 

Sihh = A 1 S^{uxU2)Sxu z Sxu A dxdt 

Shih = A -1 Sxiux^S^SxUidxdt 

Shhi =T2 A~ 1 S\(u l u 2 )S x u 3 S fJ ,u 4l dxdt 
h«\ JoJr3 

The first sum can be estimated using Strichartz and energy estimates as 
follows: 



\Sihh\ ^,^2^ 2 \\S^(uiu 2 )\\ L i L ^\\ ^2 S x u 3 SxU4 



1 ||5'm( U 1 U 2)||liL3 SUp V \\S X U 3 (t)\\ L 2\\S X U 4 (t)\\ L 2 

» *e[o,n A >^ 

S\ u i\\l^L6\\ u 2\\l2L6 SUp ||«3(*)lli?(m) 11^4(0 ll^m- 1 ) 

te[o,T] 
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The second and third sums are similar. Using the Strichartz estimates we 
obtain 

\Shlh\ < ^ A _2 ||5 , A(niU 2 )||L2 L 2||S' M U3|| L 2 LO o||S'AU4||L°°L2 

, r fj,m(\) .. . _ ... 

2^ 71 7l\\ b ^ UlU2 >h2L% \\ U3 \\u\H(m)\\ u 4:\\vZH(m-i) 

At least one of u\ or U2, say u\, must have frequency at least A. Then we 

continue with 

^ v ^ fim(X) | 

> "I — ||ui|| L oo L 2 11^21^2^6 11^311^2^) INHy-Jiffm- 1 ) 

M «A A5m(/i) 

< ^ ^ /x 2 m(A) A^ _ 

By (54) the first fraction is less than one therefore the summation with 
respect to A and /j, is straightforward. □ 

Next we consider the wave nonlinearity. If m is as in (53) then the linear 
wave equation is well-posed in H(m), and can easily define the corresponding 
spaces U^rH(m), V^H(m) respectively DUy V H(X~ 1 m). The next result 
asserts that in effect the contribution of the wave nonlinearity is one half of 
a derivative better than the solution to the Schrodinger equation. Here we 
impose an additional condition on m, namely 

, , m(A) X 13 

63 > -g, A > n 

which guarantees that the H(m) norm is at least as strong as the H@ norm. 
Lemma 26. a) Let j3 > \ and m satisfying (54) and (63). Then 

+ \\u\\ 2 ulH0\\ A \\u^H{m)) 

b) For P > | we have 
(65) \\ P (uV A v)\\ Du ^ T;H _ h) < A T%\\ ulL2 \\v\\ ulH , 

Proof, a) By duality we have two estimates to prove. The first is 

i-T i 



(66) / / uVuBdxdt 

Jo Jm. 3 

with divergence free B. The second is 
uuABdxdt 



< A T \\u\\ ulH{m) IMIt/2 H 4B\\ v , H{xh/r 



(67) 



JO JR 1 



< A T 5 || U || a 2 H(m) ||n||^P||^ Hl ||i?|| y ^ (A i /m) 

+ T*||«||^||«||^||A||^ H(m) ||B||^ H(Ai/m) 
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Consider (66). Since B is divergence free, it follows that the gradient 
can be placed either on u or on u. Hence using a simultaneous trilinear 
Littlewood-Paley decomposition of the three factors we reduce the problem 
to estimating the following two terms: 

S hh i = E E I 7 £a>> Vu,B)|, S hlh = EE Vu,B)| 

A ju<A A ^<A 



We use (60) and Lemma 8 to estimate the first term: 

A e m(//) 

A «<A 1 



The bound (63) insures that the summation is straightforward if e is chosen 
sufficiently small. For the second term we use (61) and Lemma 8: 

stah <a t 5 e E ^\\ u \\u%H( m) h\\u2 AH 4m V 2 H{ ^ /m) 

which is again summable if e is sufficiently small. 

Next we turn our attention to (67). Using a Littlewood-Paley decompo- 
sition for all terms it suffices to consider factors of type 



jT 

Ai,A2,A3,A4 



/ / S Xl uSx 2 uSx 3 ASx 4 Bdxdt 
Jo JR 3 



and prove that for some S > they satisfy the bound 

(68) Ca 2 ,a 3 ,a 4 £a Ar V W ■ RHS({67)) 

We begin with a weaker bound which follows directly from Strichartz esti- 
mates, namely 

T — ~ ~ — ~ — ~ JVii n2 

I \ 1 ,\ 2 MM iCA T3\ 1 4 A 2 4 A 3 3 A 4 ||li||^2 H ,3 ||^||^2^1 ||-B|| iooL 2 

Arguing as in Lemma 8, this allows us to replace (68) with 

(69) JaU.A3.A4 A7 V • modified RHS((67)) 



where we have replaced VyyH(Xz /m) space with the similar U 2 space in the 
right hand side of (67). 

Due to (13) both wave factors are I 2 summable with respect to unit spatial 
cubes therefore it is enough to estimate the above integral on a unit cube 
Q. Also we must have A4 < max{Ai, A2, A3}. Hence we consider two cases. 

If max{Ai, A2, A3} = Ai (or A2) then we estimate 

Jai,A2,a 3 ,a 4 ^ T ,f II'S'a 1 ^||l4l3||'S'a2^||l2 LO o||S'a3-4||l3l6||S'a4-B||l°°l2 

1 

1 x -i+e A£rri(A4) „ „ |U|| ||R|| 

<ieA x A 2 A 3 — i ll«lll/2H(m)ll u lll72H/'ll A llv»Hi|l- ll r/2 H( V\^ 
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By (63) the fraction above is less than one therefore for small enough e the 
bound (69) follows. 

The second case is when max{Ai, A2, A3} = A3. Then we estimate 

^Ai,A2,A3,A4 

< T12 ||S'a 1 u|| L 4 L 6||S'a 2 u||l2l-»I|5'a3-4||l6 L 3||S'a 4 -B|| LO o L 2 
1 

s rpJ- > \+e~f3 ±+ e -/3 -± A|m(A 4 ) 

Afm(A 3 ) 

and (69) again follows. 

b) By duality we have two estimates to prove. The first one is trilinear, 

(70) f [ aVrBd.nll 

Jo Jm? 



< A T 5 ||u||^ L2 |M|^|B|| y ^i 



with a divergence free B. The second one is quadrilinear, 
(71) 



uuABdxdt 



<A T 5 \\u\\ ulHl3 \\vWu*!? P||^i||B|| y ^i 



Since i? is divergence free, in (70) we can place the gradient on either u 
or v. The argument is similar to the one in part (a), using (60) and (61), as 
well as Lemma 8 in order to substitute the V 2 norm by the U 2 norm. The 
new restriction j3 > | arises in the case when the low frequency is on B. 
Indeed, if /j, <C A then by (60) we have 

I^A,>, Vw, B)\ Za^X 1 -^-* min{l, ^}\\u\\ V 2 aL2 \\v\\u%h4B\\ v ^ 

The worst case is when /j, = A2, when the coefficient above is T s X^~~^ +e . 

The estimate (71) is also proved as in part (a). Indeed, by Corollary 7 
we can substitute the V 2 space by U 2 at the expense of losing e derivatives. 
Because of the finite speed of propagation for the wave equation, see (13), 
we can reduce the problem to the case when A, B are supported in a unit 
cube Q. There we can use the Strichartz estimates for the wave equation, 
respectively the local Strichartz estimates for the Schrodinger equation. 

□ 

The next step in the proof of the theorem is to establish an a-priori 
H 1 estimate for H 2 solutions. This is obtained in terms of the conserved 
quantities in our problem, namely E and Q. 

Proposition 27. Let (u, A) be an H 2 solution for (3) in some time interval 
[0,T ] with T < 1. Then 

\\ u \\u\{o,T -m) + \\ A \\uUo,T -m) < c{E,Q) 



Proof. We use a bootstrap argument. Since u,Ae L°°H 2 , we can easily 
estimate the wave nonlinearity and obtain DA € L°°H l . This implies that 
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A £ UyyH 1 , and, in addition, that the function 

T Plli^(o,T;fri) 

is continuous and satisfies 

J™ll A lll^(0,T;Hi) = W)\\m + ||A t (0)|| L2 

A similar argument applies in the case of the Schrodinger equation. 
By (9) we estimate in the Schrodinger equation 

\\ u \\u\(o,T-m) % \\u(Q)\\m + ||*«t - & A u\\ DU 2 {0jT . H iy 
Then by (62) we obtain 

\\ u \\u%(o,T-,m) < Ca(\\u(0)\\ h i +r*||u||^ (0)T . Hl) ) 
Similarly we can use (64) to obtain a bound for the wave equation 
\\Mub(0,T;Hi) < PWIIiJi + Ut (0) \\ L 2 + T 5 C\ ||u||^ (0)T . H i ) 

We multiply the first equation by c\ = (C^) -1 and add to the second 
equation to obtain 

Pll Ul (0,T;ffi) + CAh\\u%(0,T;Hl) < \H0)\\m + M(0)||ffi + M t (0)|| L2 

+ T 5 C\{1 + 11^11^2(0^.^1)) 

We make the bootstrap assumption 
c\\\u\\ul{o,T-m) + \\Mu^(o,T-,m) < 2 + ||«(0)|| H i + \\A(0)\\ H i + p t (0)|| L 2 
Then the previous bound implies that 

ca\H\u\{q,t-,hi)+ UWu^T.m) < \H0)\\m + l|A(0)|| H i + p t (0)|| i2 
+ T 5 C{E,Q){l + \\u\\l l[QT . Hl) ) 
This shows that for T < T (E, Q) we have 

(72) ||«||l^(o l T;Hi) + ll^llt^(0,T;Hi) < 1 + ll«(0) \\m + P(0) ||ffi + \\A t {0) \\ L 2 

improving our bootstrap assumption. 

Hence a continuity argument shows that (72) holds without any bootstrap 
assumption. The conclusion of the proposition follows by summing up with 
respect to Tq(E,Q) time intervals. 

□ 

The next step is to establish an a-priori H 2 bound with constants which 
depend only on the H 1 size of the data. 

Proposition 28. Let (u, A) be an H 2 solution for (3) in some time interval 
[0,T ] with T < 1. Then 

\\ u \\u\H* + < c(E,Q)(\\u \\ H 2 + \\A \\ H 2 + Pi||tfi) 

Proof. The argument is similar to the one above. □ 
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Given the local well-posedness result in H 2 proved in earlier work [9] , we 
can iterate the argument and conclude that the H 2 solutions are global. 
Finally, our last apriori estimate is in intermediate spaces: 

Proposition 29. Letra be a weight which satisfies (54) and (63). Let (u, A) 
be an H 2 solution for (3) in the time interval [0, 1] . Then 

\\ u \\u\H(m) + Pll[/2,ff(m) < c(E,Q)(\\uo\\ H (m) + IIA)Htf(m) + 11^1 Ili^A" 1 ™)) 



Proof. The argument is again similar to the one above. □ 

In order to obtain H 1 solutions and to study the dependence of the so- 
lutions on the initial data we need to obtain estimates for differences of 
solutions. Given a solution (u, A) to (3) we consider the corresponding lin- 
earized problem 

{iv t - A A v = 2iBVu + 2ABu + <f>v + A -1 (3fa«;)u 
UB = P{vVau + vSJ A v + B\u\ 2 ) 
Our main estimate for the linearized problem is 

Proposition 30. Let (u, A) be an H 2 solution for (3) in the time interval 
[0,1]. Then the linearized problem (73) is well-posed in L 2 x Hz x H~2 
uniformly with respect to (u, A) in a bounded set in the energy space, 

(74) \\v\\ V 2 aL2 + \\B\\ uhH i < c(E,Q)(\\v \\ L 2 + \\B \\ H i + \\Bi\\ H _i) 



Proof. The conclusion follows iteratively in short time intervals provided 
that we obtain appropriate estimates for the terms on the right: 

\\2iBVu+2ABu+ ( / ) v+A-\Muv)u\\ DU 2 L2 <T 5 c(E,Q)(\\v^ i) 



respectively 



\\P(vV A u + uV A v + iB\u\ 2 )\\ Du ^ H _r 



<T s c(E,Q)(\\v\\ u2aL2 + \\B\\ u2wH i) 



These in turn follow by duality from the trilinear and quadrilinear bounds 



(75) 



BVu\U2dxdt 



JR 3 



<a T 5 \\B\ 



2 H h\\ U l\\u*Hi\\U2\\v*L* 



(76) 



(77) 



(78) 



ABu\U2dxdt 



o Jr 3 



nA 1 (uiU2)usU4 1 dxdt 

nAT 1 (uiU2)uzu^dxdt 
_i 3 



< A T 5 \\ui 11^2 H A\ u 2\\u\m \\uzWu\l* \\ u A\v%L2 



< A T S 1 1 Ui 1 1 V 2 a H i 1 1 U 2 1 1 (72 L 2 1 1 U 3 1 1 jj a H i 1 1 U 4 1 1 y2 L 2 



MAXWELL-SCHRODINGER 



27 



(79) 
(80) 
(81) 
(82) 



rT 
10 



BVuiU2<ixdt 



BVuiU2<ixdt 



<a T d \\B\ 



<a T d \\B\ 



v1h i\\Mu%h4Mu%LZ 



ABu\U2<lxdt 



JR 3 
T i 



ABu\Uidxdt 



JR S 



<a T d \\A\ 



<aT' 



U 2 H 1 



\B\ 



The quadrilinear mixed bounds (76), (81), (82) follow trivially from the 
Strichartz estimates. For (76) for instance we have 

rT, 



ABu\U2dxdt 



JR 3 



< TT2\\A\ 



L 3 L» 



\B\ 



L 4 \\ u i\\l 3 L s \\ u 2\\l°°L 2 



\B\ 



We note that there is significant room for improvement in this computation 
by localizing first to the unit spatial scale and then using the local Strichartz 
estimates for the Schrodinger equation. 

The quadrilinear Schrodinger bound (77) corresponds to the particular 
choice m(A) = 1 and j3 = 1 > \ in (62). For (78) we can write 



nA 1 \u\U2)uj,u^dxdt 
_;3 



^ ll«l«2| 



„ 3 MaUd „ 3 



~ T ' A \\ u i\\l3lA\ u 2\\l™lA\ u 3\\l3lA\ u a\\l™l 2 

< A T , 3||«l|| C7 2 jff i||u2|| C7 2 L 2||u3||[/2Hi||u4||y2 L 2 

Finally, the bounds (79) and (80) are identical since div B = and 
correspond to (70). (75) is essentially the same estimate. 

□ 

Proof of Theorem 1, conclusion: 

By Proposition 30 we can obtain a weak Lipschitz dependence result for 
H 2 solutions (ui,Ai) and (^2,^2) to (3), 



(83) 



IK - u 2 \\ l °°l 2 + \\ A i - A 2\\ tt2 1 < c(E 1 ,Q 1 ,E 2 ,Q 2 ) 

u w 

(ll(«i " «2)(0)|| La + ll^i - A 2 m\\ Hh + ||(Ai - A2)M\\ H - h ) 



We use this in order to construct solutions to (3) for H 1 initial data. 
Given (uq, Aq, Ai ) £ H 1 x H 1 x L 2 we consider a sequence of H 2 initial data 

(uq,Aq,A™) - («o, in H l xH l x L 2 
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The sequence (uq,Aq,A™) is compact in H 1 x H 1 x L 2 , therefore we can 
bound them uniformly in a stronger norm, 

Wi u O,Aoi A l)\\H(m)xH( m )xH(X- 1 m) < M 

where m(A) > A satisfies (54) and (63) and in addition 

lim A _1 m(A) = oo 

A^oo 

By Proposition 29 we obtain a uniform bound 

\W n \\L^H(m) + \\ An \\u^H(m) ^ M 

On the other hand, (83) shows that the solutions (u n , A n ) have a limit in a 
weaker topology, 

{u n , A n ) -► (it, A) in L°°L 2 x U^m 

Combining the two bounds above we obtain strong convergence in H 1 , 

(u n , A n ) (u, A) in L^H 1 x U^H 1 

In addition, -u will also satisfy the same Strichartz estimates as u n . Passing 
to the limit in the equation (3) we easily see that (u, A) is a solution. Due to 
the weak Lipschitz dependence it is also the unique uniform limit of strong 
solutions. Due to the Strichartz estimates we can bound the nonlinear term 
4>u in the Schrodinger equation as in Lemma 62. Then it also follows that 
u £ U\H l . 

The weak Lipschitz dependence (83) carries over to H l solutions, as well 
as the bounds in Propositions 27,29. Then the same argument as above 
gives the continuous dependence on the initial data. 

6. Wave packets for Schrodinger operators with rough 

SYMBOLS 

An essential part of this article is devoted to understanding the properties 
of the (95) flow at frequency A on A -1 time intervals. As it turns out, for 
many estimates the parameter A can be factored out by rescaling. This is 
why in this section we consider a more general equation of the form 

(84) iut - Au + a w (t,x,D)u = 0, u(0) = u 

which we study on a unit time scale. Here a is a real symbol which is roughly 
smooth on the unit scale. 

For such a problem one seeks to obtain a wave packet parametrix, i.e. to 
write solutions as almost orthogonal superpositions of wave packets, where 
the wave packets are localized both in space and in frequency on the unit 
scale. The simplest setup is to assume uniform bounds on a of the form 

\d^dla(t,x,0\<c a/ 3, \a\ + \/3\>k 

An analysis of this type has been carried out in [10], [11], [7]. If k = 2 
then one obtains a wave packet parametrix where the packets travel along 
the Hamilton flow. If k = 1 the geometry simplifies, and the Hamilton flow 
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stays close to the flow for a = 0; however, a still affects a time modulation 
factor arising in the solutions. Finally if k = then the a w (t,x,D) term is 
purely perturbative. 

For the operators arising in the present paper the above uniform bounds 
on a are too strong, and need to be replaced by integral bounds of the form 

f 1 \d^a(t,x\e)\ <c a/ 3, \a\ + \/3\>k 
Jo 

where t — > (x*, £*) is the associated Hamilton flow. The case k = 2 has been 
considered in [8]; as proved there, the Hamilton flow is bilipschitz and a 
wave packet parametrix can be constructed. The case k = was considered 
in [1]; then the term a w (t, x, D) is perturbative, and one may use the a = 
Hamilton flow in the above condition. 

In the present article we need to deal with the case k = 1. This corre- 
sponds to a Hamilton flow which is close to the a = flow. However, the 
term a w (t,x,D) is nonperturbative, and contributes a time modulation fac- 
tor along each packet. Given these considerations, we consider the following 
assumption on the symbol a: 

(85) sup/ \d%dP a {t,x + 2tt,Q\dt<ec a/3 \a\ + \/3\>l 

x,i Jo 

Let (xo,£o) £ ~^ 2n . To describe functions which are localized in the phase 
space on the unit scale near (a?o,£o) we use the norm: 

H$» := {f:(D- f) N f G L 2 , (x - x°) N f € L 2 } 

We work with the lattice Z™ both in the physical and Fourier space. We 
consider a partition of unity in the physical space, 

^2 4>x = l 4>x {x) = 4>{X - Xq) 

where is a smooth bump function with compact support. We use a similar 
partition of unity on the Fourier side: 

= i, = v(£ - Co) 

An arbitrary function u admits an almost orthogonal decomposition 
so that 

(86) Yl \\ u ^of H N,N < \\u\\ 2 L2 

We remark that a continuous analog of the above discrete decomposition 
can be obtained using the Bargman transform. 
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We first establish that the Hamilton flow is close to the Hamilton flow 
with a = 0: 

Lemma 31. Assume that (85) holds with a small enough e. Then for each 
(x°,£°) G R 2n and t £ [0, 1] we have 

\x t -(x°+2te)\ + \e-e\<e 

The proof is straightforward and is left for the reader; it is also essentially 
contained in [1]. This allows us to apply the main result in [8]: 

Proposition 32. Assume that (85) holds with a small enough e. Then 
for each N > the solution of the homogeneous problem (84) satisfies the 
following localization estimate: 



(87) \\u(t)\\ H 



N,N <j v | " U | 1 1 ■ 



We denote the evolution operator for (84) by S(t, s). If the initial data is 
uo = 8 X then it has a decomposition of the form 

(88) n = n ?«(°)' Ko(0)lln"-? £1 

By (32), at time 1 the corresponding solutions n^ are concentrated close 
to x + 2t£o, therefore they are spatially separated. Hence we obtain the 
following pointwise decay: 

Corollary 33. The kernel K '(1,0) of 5(1,0) satisfies 

\K(l,x,0,y)\<l 
The solution of the homogeneous equations (84) satisfies 

115(1, 0)u ||loo < ||n || L i 

If in addition the initial data is localized at some frequency A, say uq = 
S\5 X then the decomposition in (88) is restricted to the range |£o| ~ A. Then 
the corresponding solutions travel with speed O(A), and we can obtain better 
pointwise decay away from the propagation region: 

Corollary 34. The kernel K x (l, 0) of S{1,0)S X satisfies 

(89) \K(l,x,0,y)\<(X + \x-y\)- N , \x - y\ ^ X 
The kernel K\(t, s) of S(t, s)S\ satisfies 

(90) \K(t,x,s,y)\ < X~ N , \x-y\^X,\t-s\<^l 

The next result concerns localized energy estimates. 

Corollary 35. For each ball B r of radius r > 1 the solution u to (84) 
satisfies 

(91) \\S(t,0)S x u \\ L2{Br) < A-5r5||uo|| L2 
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Proof. We consider the wave packet decomposition for u = S(t,0)S\u(0), 

|&|«A 

Let Xr be a cutoff corresponding to B r . Since r > 1 it follows that the 
functions Xr-^x ,Co are almost orthogonal, therefore it suffices to prove the 
estimate for a single packet. But a single packet is concentrated near a 
tube of spatial size 1 which travels with speed O(X). This tube intersects 
the cylinder [0, 1] x B r over a time interval of length X~ 1 r. The conclusion 
easily follows. □ 

To obtain any results below the unit spatial scale we slightly strengthen 
the condition (85) by adding a weaker pointwise bound 

(92) \^a(t,x,0\<c a/3 (0K Va,/3 

This will guarantee that on a unit spatial scale the flow in (84) is a small 
perturbation of the flat Schrodinger flow. Then we have: 

Proposition 36. Assume that the conditions (85) and (92) hold. Then 

(i) For any r > the solution u to (84) satisfies the localized energy 
estimates 

(93) ||S(t,0)S A «(0)|| La(Br) < A-MlMtf 

(ii) For each y, z with \y — z\ pa A we have the square function bound 



(94) 



Js x S(t,s)S x (f(8)S y ) 



^ 11/ IL- 
LS 



Proof. To prove this result it is convenient to replace the I? initial data 
space by weighted L 2 spaces. 

Definition 37. A weight m : M? n — > M + is admissible if 

\m(x, g)/m{y, r))\ < (1 + \x - y\ + |£ - n\) N 
for some real N. 

Correspondingly we define a weighted I? space 

W U \\h(m) = Yl IM S 0,CoK ,&& JV ' JV 

Given a weight tuq at time we evolve it in time by 

m t (x + 2t£, ^) = mo(x, 
As a consequence of Proposition 32 we obtain 
Lemma 38. Assume that (85) holds with a small enough e. Then 

\\S(t, s)||L2( ms )-L2( mt ) < 1 
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Next we consider truncated solutions on a unit spatial scale. Given a unit 
ball B and an associated cutoff function \ we have the following weighted 
local energy estimates: 

Lemma 39. For any solution u to (84) we have 

Proof. We begin again with a wave packet decomposition of u, 

u = u *o,Ho 

(x ,Co)GZ 2 " 

The functions X u x ,£ are almost orthogonal in L 2 therefore the bound for 
Xu follows. On the other hand we have 

(ify - A)xu = -2VxVu - Axu + xa w (t, x, D)u 

The first two terms are estimated using the bound for x u - F° r the last one 
we note that, by (92), the operator a w preserves the H N > N spaces, 

\\a w (t,x,D)u\\ H N-i,N-i < (^o) 2 ||m||^]v,jv 
Hence we can use orthogonality again. □ 

Now we can conclude the proof of the Proposition. For the local energy 
estimate (91) we first truncate u to a unit scale. By the above lemma with 
m = (1 + A -3 (£) 3 )(l + A 3 (£)~ 3 ) we obtain 

^llxi«IUa(m') + *~*ll(*fls - A)xm|| L 2 (m ,) < ||u || L 2 (m) 
where m! = (1 + A~ 2 (£) 2 )(l + A 2 (£)~ 2 ). It remains to show that 
X^r-^\\x r u\\ L 2 < A5||xiu|| L 2 (m ,) + \~?\\{id t - A)xi«|U 2 (m') 

Then we can localize the right hand side to the A -1 time scale. On the 
A -1 time scale we can use the Duhamel formula to further reduce the prob- 
lem to a corresponding estimate for solutions to the homogeneous constant 
coefficient Schrodinger equation, namely: 

X^r~^\\x r e~ ltA uo\\ L 2 < ||u || £ 2 (m /) 

After a dyadic frequency decomposition this becomes 

X^r-^\\x r e~ ttA SxUo\\ L 2 < \\u \\ L 2 



which is exactly the local energy estimate for the homogeneous constant 
coefficient Schrodinger equation. 

Consider now the square function bound. For \t — s\ <C 1 we can use the 
kernel bound (90). Hence without any restriction in generality we assume 
that t E I, s E J where /, J are intervals of size O(l) with O(l) separation. 
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Choose to the center of the interval between I and J. We factor the estimate 
in two and prove the dual estimates 



S(t ,s)S\(f(s)6y] 

i 



;$a h/ik- 

L 2 (m) 



respectively 

\\(S x S(t,t )u)(z)\\ L2 < A-*|M| L 2 (m) 
where the flow invariant weight m is given by 

m(x,0 = (1 + A" 1 ]? A (x - y)\) K (l + A" 1 ^ A (x - z)\)~ K 

with K large enough. These are dual bounds therefore it suffices to prove 
the second one. 

If X is a smooth approximation of the characteristic function of B(z, 1), 
then by (a slight modification of) Lemma 39 it remains to show that v = 
xS\u satisfies 

\^\\v(t,x)\\ L 2 {J) < X^\\v\\ L 2 L 2 (m . ml) + X~^\\(id t - A)v\\ L 2 L 2 (m . ml) 

where the additional weight m' = (1 + A -2 (£) 2 )(l + A 2 (£) -2 ) can be added 
due to the localization to frequency A. 

This estimate can be localized to the A -1 timescale. In addition, since v 
has support in B(z, 1) we can freeze x = z in m and replace m by 

m(i) = (l + \- l \^{y-z)\) K 

Assuming y — z = 0{\)e\ we get 

m(0 = (l + |^ / |+A- 1 |6l) K 

Then the x' variable can be factored out and we are left with a bound for 
the one dimensional Schrodinger equation, 

||e- itA «o(-,0)|| £2(J) < \~^\\vQ\\ L 2 {ml) 
But this is exactly the one dimensional local energy estimate. 

□ 

7. The short time structure 

In this section we consider a paradifferential approximation to the mag- 
netic Schrodinger equation (18). Precisely given a dyadic frequency A we 
consider the evolution 

(95) iu t - An + iiA^VSx + 5 A ^ <v ^V)n = 0, «(0) = u 
where 

A <Vx = S <Vx A 

The multiplier S\ is added here for convenience. It guarantees that waves 
at frequencies away from A evolve according to the constant coefficient 
Schrodinger flow, thereby strictly confining the interesting part of the evo- 
lution to frequency A. In addition, the above expression is written in a 
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selfadjoint form, which guarantees that the corresponding evolution opera- 
tors S(t, s) are L 2 isometries. 

Later we will prove that on the time scale A -1 the evolution of the A 
dyadic piece of a solution u to (18) is well approximated by the evolution in 

(95) . Here we establish dispersive type estimates for (95). Our main result 
concerning the flow in (95) is as follows: 

Proposition 40. Let u\ be the solution to (95) with initial data uq,\ lo- 
calized at frequency A. Then for any interval I of size less than A^ 1 the 
following estimates hold: 

(i) the full Strichartz estimates 

(96) \\u\\\lv(i,li) S IK,a||l2 

(ii) the square function estimate 

(97) IKHz,4 (L 2 (J)) <a A~3||u 0iA ||l2 

(ii) the localized energy estimate: for any ball B r of radius r > we have 

(98) \\u\\\L*(IxB r ) %A r~5\~~5\\u 0: x\\L2 

The equation (95) is L 2 well-posed, therefore we can define the spaces 

As a consequence of the above proposition we obtain 

Corollary 41. Assume that A G U^H 1 . Then for any interval I of length 
< A -1 and any function u\ localized at frequency A the following embeddings 
hold: 

(99) \\u\\\lp(i,li) Si IKHc/2 L 2 

(!00) IK||L4 (i 2 (/)) < A A~* ^i? 

(!01) ||«A||L2(7xB r ) %A r5A~2" ||« A ||^a £2 

The first step in the proof of Proposition 40 is to establish a wave packet 
parametrix and a wave packet decomposition for solutions to (95) on the 
A -1 time scale. This is done by rescaling starting from the results in the 
previous section. 

We begin by writing in the Weyl calculus 

i(A <y rf7S x + S X A <VX V) = a w (t, x, D) 
Then the symbol a(t, x, £) can be expressed as a principal term plus an error, 

a(t, x, £) = ao(t, x, £) + a r (t, x, £) 
where the principal part ao is given by 

ao(t, x, £) = -2iA <y rx(t, x) ■ £§\(0 
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By Sobolev embeddings we have the following pointwise bound for the 
truncated magnetic potential: 

This yields 

(102) \d^a (t,x,0\ <A CapX^-W 

In addition, by the Weyl calculus it follows that a r is also localized at fre- 
quency A and satisfies 

(103) \d^a r (t,x,0\ <A CopX^-M 

This brings us to our main integral bound for the symbol a, namely 

Lemma 42. Assume that A £ U^H 1 with div A = 0. Then the above 
symbol a satisfies 

(104) ^r^^ + ^m^i^j!^ 

x£ Jo { c a p(TX)zX 2 \a\ > 1. 

Proof. The bound for a r follows directly from (103), therefore it remains to 
consider do- Furthermore, it suffices to consider the case \a\ = 0, 1, (3 = 0. 
Then we need to prove the bounds 
rT 

\A <v - x (t,x + 2tC)\dt<Th-hnX\\VA\\ U 2 L2 , |£| « A 



JO 



■/ 



/o 

respectively 

|V^ <vx (t,x + 2^)|di<T5||V^|| [/ 2 L2 , |£| ~ A 

io 

These in turn follow by dyadic summation from 

(105) sup/ \(S fl B)(t,x + 2tO\dt<T l 2 f iX-^\\B\\ U 2 L2 , |£| ~ A 

The line y = x+2t£ moves through a unit spatial cube in a time A -1 . But, 
due to the finite speed of propagation for the wave equation, see (13), the 
contributions from different spatial unit cubes are square summable. Hence 
by Cauchy-Schwartz it suffices to prove the above bound for T < A -1 . By 
(12), for T < A" 1 we have 

ll-^llt/^O^L 2 ) ~ ||-B||t/2(0,T;L 2 )> 

therefore it is enough to prove: 

[ T |(5 M B)(t, x + 2t0\dt < (AT)^A- 1 ||5||^ £2 
Jo 
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It suffices to prove the bound when B is an U 2 L 2 atom. By Cauchy- Schwartz 
it suffices to consider a single step, which corresponds to a time independent 
B. Then the last bound can be rewritten in the form 

|(^B)(z)|ds<4L|3||B|Ua 

' L 

where L is an arbitrary line segment in M 3 . We can set /i = 1 by rescaling. 
In coordinates x = (x±,x') suppose L is contained in {x' = 0}. Then we use 
Sobolev embeddings in x' and Cauchy- Schwartz with respect to x\. □ 

Next we consider the rescaling that preserves the flat Schrodinger flow 
and takes the time scale A -1 to 1, namely 

( x t 
v x (x,t)=u^ 7 =,- K 

If u solves (95) then for v\ we obtain the following equation: 

(106) ivt -Av + \- l a w D^j v = 

However, this is not sufficient, we need to repeat the same procedure for 
shorter time scales. Precisely, for each A < a < A 2 we can rescale the 
time scale to the unit scale by setting 

( x t 
v(x, t) = u I — — , — 



v 7 ^ M 

Then for v we obtain the equation 

(107) iv t -Av + a™ (t, x, D)) v = 0, a^t, x, = fi^a™ ^^v 7 ^) 

Rescaling the bounds (102), (103) and (104) it follows that this rescaled 
equation belongs to the class studied in the previous section: 

Lemma 43. For e _1 A < fi < A 2 and e small enough the symbol a M satisfies 
(85) and (92) on the time interval [0, 1] . 

This allows us to apply the results in the previous section to the evolu- 
tion (95). Rescaling the result in Corollary 33 we obtain short time pointwise 
bounds for the solution to (95): 

Lemma 44. The solution of (95) has the pointwise decay 

(108) IK*)IU<» < \t-8\-*\\u(s)\\ L i, \t-s\Ze\- 1 

Proof. W.a.r.g we can take s = 0. If A -2 < t < eA -1 then this follows 
directly from Corollary 33 applied to the equation (107) with u^ 1 = t. 

The case t < A~ 2 needs to be considered separately. For such t we split 
the evolution in two parts, 

S(t,0) = S(t,0)§ x + S(t,0)(l-§ x ) 



MAXWELL-SCHRODINGER 37 

The second part evolves according to the constant coefficient Schrodinger 
flow, hence it is easy to estimate. For the first part we use the rescaled 
parametrix in the previous section corresponding to n = A -2 . The solution 

S(t, 0)S\5 X consists of a single packet on the A -1 spatial scale which does 
not move up to time A -2 . Hence we obtain 

(109) \S(t,0)§ x S x \<X s (l + \\x-y\)- N , \t\<\- 2 

which concludes the proof. □ 

By [4], the Strichartz estimates in (96) are a direct consequence of (108). 
We continue with a decay bound away from the propagation region: 

Lemma 45. // \t — s\ < eA _1 then the kernel of S(t, s)S\ satisfies 

(110) \K(t,x,s,y)\<X 3 (l + X\x-y\+X 2 \t-s\)- N 
whenever 

\t — s\ + A -2 > A _1 |x - y\ or A _1 |x - y\ + A 2 > \t - s\ 

Proof. W.a.r.g we can take s = 0. If \t — s\ < X" 2 then we use (109). On 
the other hand if A -2 < \t — s\ < eX^ 1 and then we rescale (89) applied to 
(107) with y, = t~ l . 

Since the input is localized at frequency A, it follows that waves need 
exactly a time « A _1 |x — y\+ A -2 to travel from x to y. 

□ 

Next we consider pointwise square function bounds: 

Lemma 46. The evolution S(t, s) associated to (95) has the pointwise square 
function decay 

(111) 

Proof. If \x — y\ > 1 then we can use directly (110). If \x — y\ ^1 then we 
split the integral in two parts. If \t — s\ 3> A _1 |x — y\ then we can still use 
(110). On the other hand if \t — s\ < A _1 |x — y\ then we rescale (94) applied 
to (107) with fj, = X\x -y\~ l . □ 

We continue with the proof of (97). By the TT* argument we need to 
prove the bound 

(112) 

For this we use Stein's complex interpolation theorem. Define the holomor- 
phic family of operators 

T z f{t) = J z(t- sy-'S^t, s)S x f(s)ds 



S(t,s)S x (f(s)5 y )ds(x) 



< \r. - 



z,?m 



y\ 



-1 



I/I < eA 



-i 



S(t,s)S x f(s)ds 



<A X 2 



LiL? 



4 
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Then we need to show that 

This follows by interpolation from 

(113) ||T Z /|| L2 < ||/|| i2 , Kz = 
and 

(114) II^/ILoo^^a- 1 !!/!!^, r z = 2 

For (113) we write 

S(0, t)T z f(t) =J z(t- sY^SxSiO, s)f(s)ds 
Since S(t, s) are L 2 isometries it suffices to prove that 

z(t-sy + - 1 f(s)ds\\ L ,<\\f\\ L , 



II 

which is straightforward by PlancherePs theorem since the Fourier transform 
of zt z + l is T(z + l)(r + i0)~ z which is bounded. 

On the other hand the bound (114) is equivalent to 

\\T z (fS y )(x)\\ L2 <4 X-'Wfh* 

which we can rewrite in the form 



^(t-s) 1+i -S vx (t,s)S x (f(s)S y )ds(x) 



LIU) 



Restricting t — s to the range \t — s| < A 1 \x — y| + A 2 this is a consequence 
of (111). On the other hand for larger t — s we can use directly the pointwise 
bound (110). 

The last step of the proof of Proposition 40 is the localized energy estimate 
(98). This follows directly by rescaling from (93) applied to the equation 
(106). 

8. Short range bilinear and trilinear estimates 

We first consider L 2 bilinear product estimates where one factor solves 
the wave equation and the other solves the Schrodinger equation. 

Proposition 47. Assume that A £ U^H 1 , 1 < fi < A and \I\ < A -1 . Then 
the following bilinear L 2 estimates hold: 

(115) ||S , m (5 A Ua)|| L 2 (/xB ;3 ) <A ^ll S A||[/2 L2|K||[/2 £2, 

w A, v A 

1 

(H6) I|£a^||l2( Jx]R 3) <a ^ 2 \\Bx\\ U 2 rL 2\\u IJI/ \\ u ^ ^_ L 2 

(117) ||£^a||l2<7 xR 3) < A ^"2 11^11^^2 ||«a||[72 ^.L 2 , 
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We remark that the constants in (117) are optimal, and in effect as a 
consequence of the results in the last section (117) can be extended almost 
up to time 1. On the other hand the constants in (115), (116) are not 
optimal, but this is not so important because this corresponds to the non- 
resonant case in the trilinear estimates. 

Proof. For (115) it suffices to use Bernstein's inequality and the Strichartz 
estimates, 

\\S^B x u x )\\ L 2 < ^\\B x u x \\ s 

^ ^\\B X \\ L oo L 2\\u X \\ L 2 tL 6 

<A ^\\ B x\\u^L^\\ u \\\ul ^_L 2 

A similar argument applies for (116). 

It remains to prove (117). By (12) we replace the U^L 2 space by the 
U 2 L 2 space on a short time scale. Hence we can rewrite (117) in the form 

( 118 ) II- b m m a||l 2 (/xr3) ^a ^A _ 2 II^II^^IIvaII^ ^_ L 2 

Due to the atomic structure of the U 2 spaces it suffices to prove the above 
bound in the special case when both and u x solve the corresponding 
homogeneous equations dtB^ = 0, respectively (95). 
We consider a partition of unit on the n~ l scale 



and use the localized energy estimates (101) for u x with r = \x 



^ Yj II^O B A»II£°°II&:o«a||£2 

< A A-V^IKII 2 ^ L 2 X H^o-BJi- 
< a^V 2 !!- 6 /*!!^ 2 ^!!^!!^ 2 ^ 



A,V\ 



□ 



Next we turn our attention to trilinear estimates. We begin with the 
easier case of three U 2 type spaces 

Proposition 48. a) If \I\ < A -1 , fi < A and B^, u x , v x are localized at 
frequency fx, X, respectively A then 



(119) 



/ / B^u x v x dxdt 
JiJm. 3 



min^A 2 ) 

A w A,V\ A,^\ 
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b) If \I\ = A" 1 , /i <C A and B\, u^, v\ are localized at frequency X, \i, 
respectively A then 



(120) 



/ / B\u^v\dxdt <a/^A 1 \\B\\\u*l4 v u\\uI -L*\\w\\\u* r-i? 



Proof, a) If jU < A 2 then the conclusion follows directly from (117) since 



// 

J i Jm 



B^uxvxdxdt 

i Jm.3 



, i 



< \m\B 



Vl u M\l 2 \\ v M\L°°L 2 



l 

I- linn ti t/-\ nnn / 1 ) I i- 1~\ vnrv InJin / , 

>W 

estimate 



If n > A 2 then we use (12) to replace U^L 2 by U 2 L 2 C L 2 x Lf . Then we 



/ / B^uxvxdxdt 



WB^Wlil^WuxWl^^WvxWl^.^ 



+ 
+ 



and use the square function bounds (97) for the last two factors. 

b) In the Fourier space we obtain nontrivial contributions when either all 
three time frequencies are <C A 2 or when at least two of them are > A 2 . 
More precisely, using smooth time multiplier cutoffs we can write 

/ / Bxu^vxdxdt = X{\D t \>\y32} B \ u tiV\dxdt 

/ / X{\D t \<\2/32} B \X{\D t \>\l/32}U»Vxdxdt 
/ / X{|D t |<A2/32}- B AX{|D t |<A2/32}^X{| J Dt|<A2/16}^Adxdt 

J i Jm. 3 

Since the wave equation has constant coefficients, for the first term we 
can bound the first factor in L 2 , 

IIX{|D t |>A2/32}- B A||L 2 ^ ^ _1 ||- B A||[/2 r L2 

On the other hand for the remaining product we use the energy estimate for 
vx and the L 2 L°° bound for u„. 

We argue in a similar manner for the other two terms. The bilinear 
expressions 

X{|D t |<A 2 /32}- B AX{|D t |<A 2 /32} u M' S (i(X{\D t \<\ 2 /32}B\Vx) 

can be estimated in L 2 using (116), respectively (115). Hence it remains to 
bound in L 2 the high modulation factors: 

Lemma 49. We have 

IIX{|Di|<A2/16} w A||L2 <A A -1 ||f A ||l/2 L 2 

respectively 

\\X{\D t \>X 2 /32} u lJ.\\L 2 <A ^~ 1 \\u fJ ,\\ U ^ ^_ L 2 



MAXWELL-SCHRODINGER 



41 



Proof. In the case A = both bounds are trivial, the difficulty is to accom- 
modate the unbounded term involving A. We consider the first bound only, 
as the argument for the second is similar. 

Without any restriction in generality we can take v\ to be an U 2 atom. 
The kernels of the operators X{|D t |<A 2 /i6}> respectively X{\D t \>\ 2 /32} decay 
rapidly on the A -2 time scale. Then it suffices to prove the estimate in two 
cases: 

(i) v\ is supported in a A -2 time interval (this corresponds to steps of 
length A~ 2 and shorter). Then the bound follows directly from the energy 
estimates and Holder's inequality. 

(ii) v\ solves the homogeneous equation (95) on the time interval I with 
A~ 2 < |7| < A^ 1 (this corresponds to steps of length A -2 and longer). Then 
we can use the bound (117) to estimate 

\\(id t - A)u A || L 2 (J) = P <v ^VS , a wa||l 2 (/) In A X^\\v x \\ l ^l 2 
Hence with / = [to,ii] we can write 

(id t - A)(xiv\) = iv\(t )6t=to - iv\(ti)6t= tl + fx 
where xi is the characteristic function of I and 

II/a||l2 Si In A A5||va||l^l2 
Hence working with the constant coefficient Schrodinger equation we obtain 

IIX{|D t |<A2/2}X/^A||L2 S^ _1 (IMio)llL 2 + IM*i)||l2) + A~ 2 ||/|| L 2 

^aA _1 ||va||l°°l 2 
which is exactly what we need. 

□ 



□ 



Finally we turn our attention to the case when one of the three U 2 spaces 
is replaced by a V 2 space: 

Proposition 50. a) Let \L\ < A -1 , fi < A and B^, u\, v\ localized at 
frequency fx, X, respectively A. Then 



(121) 



B^u\V\dxdt 



<A 



(A|/|)^ 
A 1 



B Au 2 lAWxWuI -lA\ v \\\vI ^l 2 

w A,\/\ A,VA 



// in addition A 2 <c /i < A then 



(122) 



// 

J I JR ;: 



B^u\V\dxdt 



<A 



Vx 



\ B /J-\\u* L 2 \\ u x\\u 2 ^L 2 \\ v \\\v 2 ^ 
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b) If \I\ = A" 1 , fi <C X and B\, u^, v\ are localized at frequency X, [i, 
respectively A then 



(123) 



/ / B\u^v\dxdt 
J I Jr 3 



fi* In A 

A w ^.x/i 1 A, VI 



Proof, a) Using the bilinear L 2 bound (117) for the product of the first two 
factors we obtain 



B^u\V\dxdt 



< 

rZA : \\&h\\u2L 2 \\ u M\U 2 ^L 2 \\ v \\\L°°L 2 

A w A,V\ 



Then (121) follows due to the trivial embedding V 2 ^-L 2 C LfL 2 .. 

On the other hand the LHS of (122) can be estimated either as above or 
as in (119). Then (122) follows from the decomposition 

We can factor out the (95) flow by pulling functions back to time along 
the flow. Then the above relation becomes 

V 2 L 2 C lncr U 2 L 2 + a' 1 L°°L 2 , a > 1 

This in turn is true due to Lemma 8. 

b) This follows from a similar argument to the one above and by using 
(116) and (120). 

□ 

9. The short time paradifferential calculus 

In this section we prove that, given a dyadic frequency A, the evolution 
of the A dyadic piece of a solution u to (18) is well approximated by the 
evolution of the paradifferential equation (95) on time intervals of size A -1 . 
We also introduce different paradifferential truncations 

(124) iu t - Au + i(A <u VS x + S x A <u \7)u = 0, u(0) = u 

and show that they all generate equivalent spaces. The spaces associated to 

(124) are denoted by U\ vX L 2 , U\ vX L 2 , respectively U\ vX L 2 . We refer to 
the above evolution as the (A <V ,X) flow. 

A special case of the above equation is when A <v is replaced by A^A- We 
refer to that as the (^4<a, A) flow. By a slight abuse of notation we denote 
the corresponding spaces by U\ x X L 2 , etc. 

Proposition 51. a) For any interval I with \I\ < A" 1 the solution u to 
(18) satisfies 

(125) I|Saw||[/2 %a ||uo|Ua 
b) In addition, for any y/X < v <C A we have 

(126) |Mlt/2 ^(/;L2) -A \W\u2 >i/iA (J;L2) 
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The bound (125) transfers easily to U 2 spaces: 
Corollary 52. For any interval I with \I\ < A" 1 we have 

Combining this with (99) we immediately obtain the bounds (20) in part 
(i) of Theorem 9: 

Corollary 53. The solution u to the homogeneous magnetic Schrddinger 
equation (18) satisfies the Strichartz estimates (20). 

The rest of the section is dedicated to the proof of Proposition 51. 

Proof, a) From the equation (18) we obtain the following equation for S\u, 

[i dt _ A + iA^VSx + iS x A <vx v) S x u = f x 

where 

fx = S x (2iA > ^Vu + A 2 u) +i[S x , A^SxVu 
Then we have 

ll-S^llt/ 2 Ai;Lt) % INHl2 + WfxWou 2 Ai;L 2 ) 

j4 , v A A, v A 

The estimate (125) follows if we establish that the inhomogeneous terms fx 
are uniformly small, 

II/aIIdC/2 _(7;L2) (A|/|) 5 sup||5 A /'u|| [7 2 ( /;L 2) 

We consider the terms in fx- For the first term by duality we need to prove 
that 

A,sf\ X' A,VJ7 K ' ' 

We take a Littlewood-Paley decomposition of the first two factors and esti- 
mate each dyadic piece. There are several cases to consider for the integrand: 

a) Sf.AVSxuSxv with \f\ < /x < A. Using (121) yields a constant 

|/|^=(A|/|)I(^|/|)^(^) 3 

which is favorable if |/| < fi~ 2 . On the other hand using (122) yields a 
constant 

which is favorable if |7| > p^ 2 . 

b) SxAVS^uSxv with jj, <C A. Then using (123) yields a constant 

//5A _2 lnA < In A 
and a power of |/| can be easily gained as in Remark 10. 




■VuSxvdxdt 
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c) S U A\7 S u uS\v with v S> A. Then we can use (123) but only on v 1 
time intervals. We obtain a constant 

i i 
A 2 In z> 

and again a power of / is gained as in Remark 10. 

For the second term in fx by duality we need to prove that 



A uSxvdxdt 



< A {\\I\) 5 \\V\\ V 2 ^ L 2 SUP ||S>||t/2 _(/;L2) 



A, •/! 



We consider the corresponding dyadic pieces 

/ / S Xl ASx 2 ASx a uSxvdxdt 
J i it 3 

< (A|/|)T2 \\Sx 1 A\\ L 6 L 3\\Sx 2 A\\ L 6 L -4Sx 3 u\\ L ^ L 2\\Sxv\\ L i L : i 

5 —I _ 2 

< A (A|7|)T5A 1 5 A 2 5 ||^||^2 ffl ||5 A 3U|| L oc L 2||«|| y 2 r2 

W A, v A 

where for u we have used the short time Strichartz estimates. The summa- 
tion with respect to Ai and A2 is trivial. So is the summation with respect 
to A3 since the integral is zero unless either A3 = A or A3 < max{Ai, A2}. 
Finally for the commutator term in fx we have the bound 

\\[Sx,A <vx ]VS x u\\ L i L 2 < \-\\\I\)\\[Sx,A <vx }VSxu\\ L oo L2 

<X-\\\I\)\\VA <vx \\ L «>\\S x u\\ LooL 2 

<A-i(A|7|)P|| L oc^||5 A n|| LO c L2 

which again suffices by duality. 

b) By the same argument as in part (a) we obtain 

W A <V>y u ~ a <» Vu Wdu 2 ai;L*) (A|/|) ,5 ||u|| [7 2 (/;L2) 

which shows that the two flows are close. Then (126) follows due to Lemma 3. 

□ 



10. Generalized wave packet decompositions and long range 

trilinear estimates 

Denote by Ci(fi, X,\I\) the best constant in the estimate 



/ / B^u x vxdxdt 
J 1 JR 3 



< Ci(fi, A, |/|)||S m ||^ L 2||ua||{72 iAiA l2||ua||{72 iAiA l2 



with B^, ux and vx localized at frequencies /j,, A, respectively A. Similarly, 
let C2(/x, A, |/|) be the best constant in the estimate 



// 

J I 



Bxu^vxdxdt 



< C 2 (n,\,\I\)\\B\\\ u * rL 2\\u li \\ u * vWwxWu* aaL 2 



with Bx, U[i and vx localized at frequencies A, /x, respectively A. 
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As a consequence of Proposition 48 we have 

Ci(//,A, A -1 ) < A A" m in{//A _ 2,l} 
A trivial summation shows that for larger time intervals we have 

(128) Ci(m, A, |/|) <a A-5 min^A^, 1}(1 + A|/|) 
We seek to iteratively improve this to 

(129) d( M , A, |/|) < A A-3 min{/xA-5, 1}(1 + A|/|)5 

for intervals I almost up to length 1. To achieve this we iteratively produce 
an increasing sequence of times T for which (129) holds for |7| < T. At the 
same time we seek to improve Cii^li-, A, |/|) in a similar manner, as well as 
extend the time for which the local energy and local Strichartz estimates 
hold. More precisely, we consider a set of properties as follows: 

(i) Paradifferential approximation of the flow: 

(!30) II-SamII^^^) <a IM|[/2(/,l2) 

(ii) Trilinear bounds: 

(131) Ci(/i,A,|/|) <4 A^min{^A-5,l}(l + A|/|)5 

(132) C 2 (/x,A,|/|)<4 ^A-HH- A|J|)s 

(iii) Local energy and local Strichartz estimates for each cube Q of size 1 
and each function u\ localized at frequency A: 

( 133 ) IM|l2(/;L2(Q)) <A A~2 ||ma|Ie/2 jA]A L2, 

( 134 ) II u aIIl2 (/;L 6(q)) <a II u aIIc/2_ a a l2. 

So far, by Propositions 40,51,48 we know that the above estimates (130)- 
(134) hold if |/| < A -1 . On the other hand, in order to prove Theorem 9 we 
need to know that (130)-(134) hold if |/| < A~ e for e arbitrarily small (see 
also Remark 10). This is accomplished in the next result. 

Proposition 54. Let < a < 1. Assume that the estimates (130)-(134) 
hold for \L\ < \- a . Then (130)-(134) hold for \I\ < A _/3 for each (3 > fa. 

Proof. We first improve the time range of the paradifferential calculus: 

Lemma 55. a) For each frequency A we have 

ex. 3 

\\ S \ u \\ul xx (I,tf) St |M| 178(7,1,2), |/| < A~2(logA) _ 3 

b) For each frequency A 1- ~ log A < v <C A we have 

IMIt/2 ^ A (7,L2) «a ||^||[/2 jA A (7,72 ) , \I\ < T(A, v) = i/A~ 1- *(]ogA) -1 
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Proof, a) We observe that (131), (132), (133) for |/| = X~ a trivially lead to 
bounds for longer time, 

(135) dfji, A, |J|) <4 A-l min^A"*, 1}(1 + + A Q |7|)§ 

(136) C 2 (/i, A, |/|) <4 /iU-^l + A|7|)3(l + A Q |/|)5 

(137) H«a||l3(j,l6(0)) <a (1 + A Q |/|)5 IIuaII^^ 

Furthermore, due to (130), we obtain the same constant in the trilinear 
estimates when we replace u x by S x u and ||«a||i/2 a l2 by ||ti|| c/^z, 2 • Also by 

the argument in Lemma 8, we can also replace one of the U 2 norms with a 
V 2 norm at the expense of an additional In A loss. 

The rest of the proof is similar to the proof of Proposition 53. For u 
solving (18) we write 

(id t - A + iA <<x VS x + iS x A <<x V^J S x u = f x 

where 

fx = S x (2iA> x Vu + A 2 u) + i[S x , A <<x ]S x Vu 
Hence it suffices to prove that 

We use duality and consider each term in f x . For the first one we need to 
show that 

T J^ A >x^uS x vdxdt < A \\Mu^h4 u \\u\(i,l2)\\ v \\vI xa lz 

After a Littlewood-Paley decomposition of the first two factors we need to 
consider the following three cases for the integrand: 

(i) S x A\7S x uS x v. Then we use (135) to obtain a constant 

logA A3(A|/|)5(A a |/|)3A- 1 = logA A§|/| 

which is satisfactory given the range for /. 

(ii) S x AVS^uS x v, n <C A. Then we use (136) to obtain a constant 

/i2A- 1 (A|/|)5(A Q |/|)U- 1 < At \I\ 

which is much better than we need. 

(iii) S u AVS u uS x v, A <C v. Then we use (136) to obtain a constant 

Afi/ _1 (i/|7|)5(^ Q |J|)§^- 1 < At \I\ 
For the second term in f x by duality we need to prove that 



ABuS x vdxdt 

I 



WMu^m WBWu^hi ll^llf/2 (/;i2) ||u|| v a a aL 2 



Due to the finite speed of propagation for the wave equation, see (13), it 
suffices to consider the case when A and B are supported in a unit cube. 
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Then we bound A and B in L°°L , u in as in (137), and v in L°°L 2 and 

use Holder's inequality with respect to time. This yields the same constant 

a 

A"2|i] as above. 

Finally we consider the commutator term in f\. This can be represented 
in the form of a rapidly convergent series of the form 

[S X , A <<x ]VS x u = £ Sl 3 (VA <<x S 2 x 3 u) 
j 

where , are operators similar to S X - 
Then by duality it suffices to prove that 



V A^x^ S\uS\vdxdt 

J./l»2 



3 a 

<A (log A) a A 2 | J| || VAHtfa^a IMI^a ||v|| y* 



For this it suffices to consider a Littlewood-Paley decomposition of A^ x and 
to apply (135) for each dyadic piece. 

b) By virtue of Lemma 3 it suffices to show that 

WA^x^Sxu- A <u VS x u\\ DU 2 a x x < A \\u\\u2 a x x 
and the similar bound for S^A^aV — SxA <u V. By duality this becomes 

A 



< A \\Muih4u\\uiJ\v\\vi XiX 



/ / SuAVSxuvdxdt 

Indeed, for |/| > X a the estimate (135) yields a constant 

logA A5(A|J|)3(A a |J|)3i/- 1 = logA zz-^f+Vl 

which leads to the restriction 

|/| < T(A, v) = {\og\)- l vX~^- 1 

We observe that this is useful only if it provides information on time intervals 
with |7| > X~ a . This leads to the condition 

v > A 1- ^ logA. 

□ 

Next we consider the (124) evolution and we construct a generalized wave 
packet structure for the flow. The frequency scale is <5£ = v and the time 
scale is T(A, v) therefore it is natural to define the spatial scale by Sx = 
vT{\, v), as in the case of the flat flow. 

We first partition the initial data. Let <p be a smooth unit bump function 
in R 3 x R 3 so that 

J2 4>{x-k,t-j) = i 

Denote 

am i c\ a,( x ~ x o ?-6 s 
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where 

(xoiCo) £ 2^ x3 = (isT(\,v)Z) 3 x (z/Z) 3 
Then consider an almost orthogonal decomposition of the initial data 

Oo,£o)GZ? x3 

We denote the corresponding solutions to (124) by u XQ £ , and we call them 
generalized wave packets. To measure their evolution we consider the family 
of operators 

Wo = K^-Co), (isT(\,is))- 1 (x-x -2tt)} 

which commute with idt — A. The following lemma shows that u XOj £ is 
concentrated in a tube 

K hio = = \x - x - 2t£>| < vT(X, v), |£ - £ | < v} 
and decays rapidly away from it. 

Lemma 56. The solutions u xo £ for the (124) flow satisfy 
(138) 

E E ll^o^oWII^^) ^ ||no||| 2 , |/| < T(X,u) 

(zo,?o)ez£>< 3 \a\<N 

Proof. At time we clearly have 

E WK^o^mh <ho\\h 

\a\<N 

therefore it suffices to prove that a single generalized wave packet satisfies 

(139) £ \\L^ u X0 ^f ulM , L2) <A E HWo*Wo(0)Hi 2 

|a|<JV \a\<N 

This follows iteratively from 

(140) llA^o^l^^/,^) ~A \\L Xo ,z v(0)\\ 2 L 2 + \\ v \\ll vX {i^) 
for which, in turn, we need the commutator bound 

(141) \\[L X0 £ ,A <v S X V]v\\ DU 2 ^ A (7, L 2) ||v||^a vX (i,l*) 

as well as the similar one for the operator S\A <U V. 
If L xo£o = ^ _1 (C-^o) then 

A<„SaV] = ^(VA^^V 

therefore by duality we need to show that 



/ / VA <u X7S x uvdxdt 
J i Jr 3 

which follows from (135). 



< a ^ a \\u%h4AuiJHvi vA 
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If L XQ £ = (i/T(A, v))~ l {x -xq- 2i£) then 

[L X0 £ ,A <V S X V] = (vT(\v))- l (A <v + t(VA <v )V)S x 
The second term is as above. For the first by duality we need to show that 

which is much weaker and follows again from (135). 

□ 

The parameter v is chosen so that the packets move away from their initial 
support by the time X~ a . For this we impose the condition 

vT{\v) < A 1- " -6 

for some e > 0. This is satisfied if we choose v of the form 

u = X 1 -^ 

in which case we have 

uT{\u) = A-T-^logA)" 1 

This leads to the choice of (3 in the proposition. Indeed, the proof of the 
proposition is concluded due to the following lemma: 

Lemma 57. Let e > 0. Choose v so that vT(\,is) < \ l ~ a - £ and v < A 1_£ . 
Then (131)-(134) hold for \I\ < T(A, v). 

Proof. By the previous lemma we can replace the {A<^\,\) flow by the 
(A <u , A) flow in (131)-(134). We begin with (131). Without any restriction 
in generality we can assume that u and v are U\ vX L? atoms. The gener- 
alized wave packet decomposition in Lemma 56 easily extends to U\ v X L 2 
atoms. Indeed, if we denote by S v (t,s) the evolution generated by (124) 
then an atom u\ of the form 

k 

can be partitioned as 

( 142 ) « = E E u *o,z 

xo |?o|~A 

where 

( 143 ) ux £ {t) = Y J hk,t k+ ,)S.{t,t k )u k X0 ^ 

k 

with 

Then, denoting 

iK,,di 2 = E E \\ L %^X^A 

k \a\<N 




A <u S\uvdxdt 
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we have the orthogonality relation 
(144) 



E E !- 



Eii Ufc iii 2 



< 



We argue in a similar manner for v. Hence it suffices to take u\ as in 
(142), (143), and similarly for v\, and prove that for |/| < T(A, v) we have 



/ / B^uxvxdxdt 
J I Jk? 



< A RHS(131) ■ L 2 



(145) 



EE lll^o^olll EE 

III v xq, Colli 



xo |^ol~A 



x o |?o|~A 



We proceed with several reductions, which will eventually lead to shorter 
time intervals. 

1. Reduction to spatial scale AT(A,z^). Heuristically, in time T(X,u) 
the frequency A waves for (124) travel by AT(A,z^). Hence we partition the 
space into cubes {Qj}j^z 3 °f si ze v). Correspondingly, we decompose 

u\ into 

= E u ^ u i = E Ux o,Ho 
iez 3 x eQj 

and similarly for v. By (138) the functions Uj decay rapidly away from an 
enlargement of Qj. Precisely, if xo £ Qj and \ j — k\ > 10 then the separation 
between the tube T£ Q ^ and Qk is 0(\j — k\XT(X, v)). Comparing this with 
the tube thickness vT(\, v)) we have 



\Mt,x)\<\- N \j-k\ 



-N 



E E in Ux o '?oii 2 ' x e Qk, \j - k\ > 10 

xo£Qj |£ |~A 



Thus in (145) it suffices to consider the output of Uj and Vj 1 for \j — ji\ < 20, 
and only within an enlargement CQj of Qj; the rest is trivially estimated 
using the above bound. Furthermore, by Cauchy-Schwartz it suffices to 
consider a fixed j and fc. By a slight abuse of notation we set k = j in the 
sequel. Then (145) is reduced to 



/ / XCQjB^UjVjdxdt 

Jl JR3 



<4 RHS(131) ■ II^H- L 2 



(146) 



E E 

K x £Qj |£o|~A 



E E l v *o£o I 
.XQ&Qj |$ |~A 



2. Reduction to small angles. Here we partition the A annulus A\ in 
frequency into small angles of size ^ with centers in 9 C § 2 , 

Ax = \J A X) 8 
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Then we divide 

oee x eQj £eA e 

and similarly for Vj. It remains to prove that 



// 

J I JR ; 



XCQjBpUjfiVj^dxdt 



<4 RHS(131) ■ L2 



1 i 

2 / \ 2 



(147) 

E E EE 

lll^o^o 

K x eQj £o£A x ,e J \xoeQj (o^A x 

3. Reduction to a spatial strip of size X £ vT(X,v). Given directions 
9 and u as above we choose a coordinate, say £i, so that both the 8 and the 
ijj sectors A$, are away from £i =0. Dividing the spatial coordinates 
x = (xi,x') we partition the space into strips Sk of thickness X £ uT(X, v) in 
the xi direction. 

Arguing as in (13), the A factor is square summable with respect to the 
strips Sk- There are about X 1 ~ e v~ 1 such strips which intersect 30Qj. Hence 

by losing a A" fi~2 factor we can use Holder's inequality to reduce the 
problem to the case when A is supported in a single spatial strip Sk- It 
remains to prove that 



// 

J I JR^ 



XSknCQjB/iUjfiVj^dxdt 



< A RHS (131) ■ \-^^\\B4 u? r 2 



1 i 

2 / \ 2 



(148) 

E E 

III U X , Colli 

EE 

K x <=Qj £o£A x ,e / \xoEQj £oeA\ tU 

4. Reduction to spatial scale X £ vT(X,v). Due to the choice of co- 
ordinates above, the packets in uj a and Vj jU) travel in directions which are 
transversal to Sk- Hence if we partition Sk into cubes Qi of size X £ uT(X, u), 
each packet will intersect only finitely many cubes. Then we partition fur- 
ther 

u j,o = ^2 u j,e,i, u j,e,i= ^2 Ux »£o 

where 

A jA i = {(x ,^);x G Qj, Co G A x>0 , {x +M£ } n Qi + 0} 

and packets u XQ £ intersecting more than one cube are arbitrarily placed in 
one of the terms. 

Arguing as in the first reduction, for | Z — Zx | 3> 1 the size of Ujfij in is 
rapidly decreasing, 

Ke,ifa*)\ 2 Z*~ N \l-h\- N l««o,&l 2 » x€Q h , |Z-Zi|»l 
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Hence the joint contribution of Uj t $ t i and Vj^^ in (148) is nontrivial only on 
the diagonal |Z — Zx | ^ 1- By Cauchy-Schwartz with respect to I it suffices 
to estimate this contribution for fixed l,h, in an enlarged cube CQ[. By a 
slight abuse of notation we set I = l±. Then we need to show that 



/ / XCQiB^Uj^iVj^idxdt 

J I JR3 



<4 RHS(131)-X-^^\\B4 V 2 L2 



(149) 

S 11^0,^0 IIP I I Yj Ko,&I 

5. Reduction to a time interval of size A e "Vr(A,i/). Each tube 
T Xo £ intersects the cube Qi in a time interval of size X £ ~ 1 uT(X,v). Hence 
it is natural to partition the time interval / into subintervals I m of length 
X £ ^ 1 uT(X, v). Correspondingly, we split Ujfi,i into 

u j,9,l = ^2 u j,e,l,m, Uj,e,l,m= ^ Ux ,€o 

m (zo,£o)eA,-,<M,™ 

where 

Aj,e,i,m = {( x o> £o); x o ^ Qj, £o e A A)fl , T£ 0)& n/ m xQi/l} 

and similarly for Vj jU ,i- Again, the size of ujjj^ in I mi x Qi is negligible if 
\m — mi| S> 1. Thus by Cauchy-Schwartz with respect to m the estimate 
(149) reduces to the case of a single interval Im, 



f f 



XCQiB^u jfi ^^ m v jyU)y i ym dxdt 



< A ##5(131^-^11^11^2 



(150) 



i. i 

2 / \ 2 



^ III^O^O III I I III^O^oll' 



^(a;ojSo)e-Aj,e,j,m / \(zo,£o)eA,>,i,m 

But this follows from the hypothesis since 

|J m | = \ £ - l vT{\v) < \- a 

and 

1 l-E 1 

|/|2A 2 ^2 = |J m | 

In the case of (132) the argument is similar but with several adjustments 
which we outline. 

1. Reduction to spatial scale max{//, A e i/}T(A, v). This smaller initial 
localization scale is possible since frequency u Schrddinger waves travel with 
speed fi, so within time T(X,v) they can spread only as far as fiT(X,v). 
Thus for the factor we have square summability on the /j,T(X,v). For 
the wave factor B\ by (13) we have square summability on the same scale, 
therefore we are allowed to localize spatially the estimate on the /xT(A, v) 
scale. If \i is small we only can take partial advantage of this due to the 
wider spread of frequency A pakets. 
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2. Reduction to small angles. This is as before. 

3. Reduction to a spatial strip of size X £ vT(X,v). The differ- 
ence here is that we have only about max{l, /iA~ £ ^ -1 } strips intersecting a 
max{/i, A e ^}T(A, v) cube, therefore we only loose a factor of 

max{l,//2A 2} 

4. Reduction to spatial scale X £ vT(X,v). This is as before. 

5. Reduction to a time interval of size X £ ^ 1 uT(X, v). As before, the 
frequency A wave packets spend a time X £ ~ 1 vT(X, v) inside a X £ vT(X, v) cube 
Q. However, the frequency fi packets spend a longer time [i~ l X £ vT(X,v) 
inside Q. Hence we can carry out first a lossless reduction down to time scale 
min{l, fj,~ 1 X e v}T(X, v). To further reduce the time scale to X E ~ 1 uT(X, v) we 
can only use the square summability for the frequency A waves, therefore 
we apply Cauchy-Schwartz and loose an additional factor of 

min{l,/i _1 A £ z/}5(A £ ~V)-i 
Finally, combining the two losses in Steps 3 and 5 we obtain a total loss 

of 

which is identical to the one in Step 3 of the proof of (131). We conclude as 
before. 

The argument is considerably simpler in the case of (133) and (134). 
There each packet intersects Q x / in a time interval which is shorter than 
X- 1 vT(X,v) < X~ a . Grouping the wave packets with respect to such time 
intervals we obtain the square summability of the outputs and reduce the 
problem to the shorter time scale X~ a . 

□ 
□ 
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